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FOREWORD 


This final report presents the results of work performed by 
personnel of the Flight Technology Group of Lockheed’s Huntsville 
Engineering Center for NASA MSFC under Contract NAS8-33807. The NASA 
Contracting Officer’s Representative and technical monitor for this 
contract was Mr. Charlie C. Dill, Jr., to whom the authors are grateful 
for his valuable assistance, direction, and contributions to the 
successful completion of this study. 
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1. INTRODUCTION 


Research was conducted to determine the feasibility of replacing the 
Solid Rocket Boosters on the existing Space Shuttle Launch Vehicle (SSLV) with 
Liquid Rocket Boosters (LRB) . As a part of the LRB selection process, a 
series of wind tunnel tests was conducted along with aero studies to determine 
the effects of different LRB configurations on the SSLV. Final results were 
tabulated into increments and added to the existing SSLV data base. 

The research conducted in this study was taken from a series of wind 
tunnel tests conducted at Marshall Space Flight Center’s 14-inch Trisonic Wind 
Tunnel. The effects on CAF, CNF, CMF, CY, CSR, CTR, CBR were investigated for 
a number of candidate LRB configurations. The aero effects due to LRB 
protuberances, ET/LRB sep. distance, and aft skirts were also gathered from 
the tests. Analysis was also conducted to investigate the base pressure and 
plume effects due to the new booster geometries 

Section 2 of this study discusses the test results found in Phase I of 
wind tunnel testing. Section 3 discusses the results in Phase II of test- 
ing, along with a comparison to Phase I tests. Section 4 gives preliminary 
LRB lateral/directional data results and trends. Section 5 discusses the 
protuberance effects and section 6 the gap/skirt effects. Section 7 discusses 
and gives results of the base pressure/plume effects study. 


1 


LOCKHEED- HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 


2. PHASE I LIQUID ROCKET BOOSTER (LRB) TEST DATA 

This section delineates the methods and results of the Phase I (TWT0707) 
wind tunnel test. Section 2.1 presents the test configurations and conditions 
used during testing, and the methodology of determining the incremental data 
is discussed in Section 2.2. Section 2.3 presents the interpolations 
performed on the LRB data. 

2.1 TEST CONFIGURATIONS AND CONDITIONS 


Testing was conducted to determine the effects of length and diameter on 
aerodynamic coefficients. The configurations tested used diameters of 12.2 
ft, 15 ft, 18 ft, and 21 ft. Three lengths were tested for each diameter 
configuration ranging from 144 ft to 190 ft. Figure 2-1 and Table 2-1 pre- 
sent the test configurations used in Phase I wind tunnel testing. Table 2-2 
shows the test matrix for Phase I. The lengths were measured from nose tip to 
base, excluding the nozzle. The LRB/ET attach points were the same for all 
LRB configurations, and the base of the nose cones was not to extend aft of 
this point. The nose tip had a radius of 1.11 ft, and nose half angles were 
all 18 degrees. All LRB’s were tested - without protuberances. 

The angle of attack for the tests ranged from -10 degrees to +10 degrees 
in even increments . The Mach numbers used were 0.6, 0.8, 0.9, 0.95, 1.05, 
1.15, 1.25 1.48, 1.96, 2.74, 3.48, and 4.96. The sideslip angle was zero. 
Nominal settings of 10 degrees/5 degrees were used for the inboard/outboard 
eleven settings. 
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Fig. 2-1 Phase I LRB Test Configurations 
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Table 2-1 PHASE I LRB CONFIGURATIONS 
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2.2 INCREMENT DEVELOPMENT METHODOLOGY 


The development of the coefficient incremental data was initiated by the 
analysis of Phase I data obtained from wind tunnel testing. The data were 
compared to the current SSLV data and evaluated for validity. The increments 
were developed from the difference between a particular set of LRB data and 
that of the DILI configuration, the latter being equivalent to the current SRB 
less the aft skirtand protuberances. The resulting increments could then be 
added directly to the SSLV data base. Increments were developed for each LRB 
configuration and test point. 

For the LRB Phase I effort, the test Mach numbers are shown in Table 2-2. 
Incremental data at Mach numbers corresponding to 0.95, 1.05, 1.15, and 3.48 
were generated for the "D2", "Da”, and *’D4" configurations by subtracting 
experimental baseline (DILI) data from data generated by linear interpolation 
using the method outlined in Fig. 2-2. 

Analysis of the incremental data at the above transonic Mach numbers has 
shown that valid increments cannot be generated using this method. However, 
at higher Mach numbers where changes in the coefficient data with changes in 
Mach number are small, this method can be used with relatively good results. 

Shown in Fig. 2-3 is a flowchart detailing the approach recommended for 
any future incremental data base development efforts. The new method differs 
from the old method only in that the interpolation for desired Mach numbers is 
not performed until after incremental data are generated. This new method- 
ology for obtaining incremental data at desired Mach numbers (between Mach 
numbers for which experimental data are available) will assure that consistent 
increments will be generated. 

Shown in Fig. 2-4 are several example plots which illustrate the differ 
ences between Phase I incremental data obtained using the previous interpo- 
lation method and the recommended future method. Figures 2-4 verifies 


5 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



Table 2-2 MATRIX OF DATA FILES USED TO GENERATE THE PHASE I 
LRB INCREMENTAL DATA 


LMSC-HEC TR F268592 


a 

-j 

u. 

1 

-J 

Ul 

:c 

a 

e 

'M 

m 

5 




«• o 

X Ur 


— a* n » 

ssss 

5 555 

H H *- I- 

9* * • 

3 3 8 8 

1 - f- H» 

EffEE 

n •* 0 « 
n n n n 

EEEE 

SSR 

SEE 

n •» m ifl 

• 0 0 0 

ssss 

n n t 
0 0 0 0 

CO 0 1 - 
0 00 

EEEE 
• • • • 

EEEE 

■ • • • 

EEEE 

• • • • 

EEE 

• • • 

r~ » o* o 
m m <5 
o o o o 

EEEE 

<m o 0 i/» 

li n rl « 

5555 

r r r ? 

• r* 0 » 

i- |. r- i«. 
CM CM CM CM 

EEEE 

• • • • 

O t n 
0 0 0 

«M CM CM 

565 

Hhh 

• ft ft 

n » « # 

'T O If) Ift 
O O O O 

EEEE 

o* O — ■ fM 

0 i- p- f- 
o o o o 

EEEE 

0 O " N 

cm n n o 

EEEE 

O - fM 
0 0 0 

566 

r H P 

r- « cm o 

O <1 to 

N N N N 

5555 
•- 1 - 1 - r. 

• • • • 

" N (1 0 
« (0 0 0 
CM <\| CM CM 
t H ^ H 

V * M * 
r- *- H H 

• • • ft 

0 « h 0 
CO 40 CO 0 
CM CM CM CM 

EEEE 

• ■ ■ ■ 

0 O -* 
« r* t» 
CM CM CM 

EEE 

• • • 

m n « v i 
•* 0 0 0 

<M <M N CM 

5555 

hl-Pl- 

• * a a 

« 1 - 0 O* 
0 0 0 0 

CM fM CM CM 
(—♦—*—»— 

F- ♦- 1- H 

• 1*1 

o -« cm n 
0000 

CM CM CM CM 

EEEE 

« • • ■ 

0 CO 0 
CO CO cn 
CM CM fM 

5 6 5 

t- r- £ 

• ft ft 

01 O ~ cm 

-T *"> ft) 0 

o o o o 

t 4 t t 
►— * — ►— 

n » n e 
i- »■* »- »* 
O O o o 

EEEE 

o ft 

CM CM CM CM 

EEEE 

000 

0 0 0 

h 0 « O 
NNNn 
cm cm cm cm 

EEEE 

• ii» 

w N O f 

noon 

CM CM CM CM 

EEEE 

• » a • 

« « 0 

n n n o 

*M CM CM PM 

EEEE 

• » • » 

9 O *t 
ri m ft 

CM CM CM 

EEE 

« • • 

cm o 0 

N (3 N N 

EEEE 
• « • • 

0 l> 0 0 

PM CM CM CM 

Em 

• a # a 

o « n n 

CM CM CM CM 
CM CM fM CM 

EEEE 

« • • • 

0 CO 0 
CM CM CM 
CM CM CM 

EEE 

• ■ • 

W 0 •* 0 

0 0 0 0 

o o o o 

EEEE 

r* 0 at o 

f» <• I' 0 

o o o o 

EEEE 

"NO* 
CM CM CM CM 

EEEE 

0 P- 0 

0 0 0 

P pH 

EEs 

•* • o» o 

Ol Ol 91 O 

EEEE 

• • a • 

« cm n ** 
o o o o 

CM CM CM CM 

EEEE 

• • * • 

SZS* 

EEEE 

• • • • 

S2: 

CM CM CM 

• • • 

- N (1 f 

f 0 

o o o o 

!* t * £ 

i- «- 

"NOt 

® » m m 

o o o o 

EEEE 

r» 9 o> o 

iiii 

0» O 
0 O CO 

565 

»—*—►— 

1 * C 9) O 
n n n 0 
o o o o 

Ht-CC 

<o r- • » 

• 09 0 0 

h P ? P 

• • • ■ 

o « n n 

0 9 0 Ok 

EEEE 

• • ■ • 

0 CO 0 

« o» 5 

556 

f- H K 

ft ft • 

n « n « 
noon 
o o o o 

EEEE 

•n « h ■ 

• 9 • 

o o S o 

EEEE 

n « o 0 

nil 

NO* 
Cl LI O 

EEE 

0 O ~ M 

n o « o 
o o o o 

EEEE 

0 > • *- • 
p* r» r- 

EEEE 
• • • • 

• O n N 

P -000 

EEEE 

ft ft ft 0 

n 0 o 

000 

EEE 

• • • 

n « h • 
N N N N 
O O O O 

EEEE 

2SJ« 

EEEE 
• • • • 

0 * O M 

« 0 «- n 

||ii 

CM O 0 
P C*" l« 

EEE 
• • • 

« N « ^ 
N N N N 
O o o o 

SEES 

• Ohn 

• • • 0 

o o o o 

EEEE 

32 = 2 

EEEE 

0 0 *» 
CO 0 ft) 

EEE 

h a » O 
o o o o 

EEEE 

iisi 

EEEE 

8 8 S 3 

EEEE 

0 9 0 

0 0 0 

EEE 

n <r n 0 
o o o o 

EEEE 

h 00 O 
» ; « o 

O O O ~ 

1? & £ H 

" N ft * 

o O o o 

H £ H H 

0 n n 
0 0 0 

EHH 

« n n » 

-* -» -4 

Q 

" n n t 

J Ml J J 

CM 

o 

" N f) * 
-1 -1 J «J 
O 

a 

553 

5 


2 w «o i- 


o 

o» 

o 


I 


6 

LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 



Fig. 2-2 Methodology used to Generate Incremental Data 
for the Phase I LRB Effort 
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Fig. 2-3 Recommended Future Methodology that Should Be 
Used to Generate Incremental Data 
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M = 1.05, a = -4 



Fig. 2-4a Comparison of Increment Development Methods ( ACNF) 
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Fig. 2-4b Comparison of Increment Development Methods (ACSR) 
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Fig. 2-4c Comparison of Increment Development Methods (Actk) 
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that the previous increment development method used to generate the Phase I 
data is not responsible for the large differences that exist between the Phase 
I and the Phase II incremental data at transonic Mach numbers. 

2.3 VERIFICATION OF INTERPOLATIONS PERFORMED ON LRB DATA 


Once raw data were available and an incremental development chosen, the 
results were interpolated versus alpha and Mach and compared to existing SSLV 
data. Section 2.3.1 gives results for angle of attack interpolations, and 
section 2.3.2 gives results for Mach number interpolations. 

2.3.1 Angle of Attack Interpolation 

Phase I results were interpolated to even angles of attack ranging from 
-8 degrees to +8 degrees in 2 degree increments from the raw wind tunnel 
data. Attached are the plots of various coefficients as a function of attack 
(Figs. 2-5 to 2-10). As the plots show, the interpolated data overlay the raw 
data extremely well for relatively linear curves (CMF and CNF). For non- 
linear curves (CAF), there are some slight differences in the curves 
representing raw and interpolated data. 

These differences are due to using a linear interpolation on a nonlinear 
curve and not to faulty interpolation. In all cases, the differences are 
quite small and are well within the accuracy range of the experimental data. 

In future efforts, higher order interpolation methods will be used in 
interpolating non-linear and/or critical data. 

2.3.2 Mach Number Interpolation 

The raw data from the wind tunnel testing were interpolated for Mach 
numbers within the LRB Phase I data base. Table 2-3 presents the LRB Phase I 
configurations represented in the data. The attached Mach numbers plots 
(Figs. 2-11 through 2-12) detail the data that were linearly interpolated from 
the raw data. The circled data points shown in the plots lie in a straight 
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TWT0707 MACH * 0.6 



ALPHA 


Fig. 2-7 CNF vs a 
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Fig. 2-8 CNF vs a 
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Table 2-3 PHASE I LRB DIMENSIONS 


(ft ) 

Length (ft ) 

D1 

D2 

D3 

D4 

12.2* 

15 

18 

21 

LI 

144* 

149 

154 

158 

12 

150 

159 

163 

167 

L3 

170 

172 

175 

177 

L4 

190 

190 

190 

- 


★ 

Currant SSLV configuration. 
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Fig. 2-11 CNF vs MACH (<* = -4) 
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line between the actual data points. From the plots it is evident that the 
data base generated was compatible with the existing operational Shuttle 
aerodynamic data. 

2.4 EVALUATION OF PLUME EFFECTS 


An analysis was conducted to investigate plume effects on Shuttle LRB 
configuration aerodynamic characteristics. The analysis included prediction 
of the correlation between base pressure and various parameters such as the 
number of nozzles, thrust, area ratio, chamber pressure and base geometry. 

The plume effects were evaluated with regard to axial force, normal force and 
pitching moment. 

The method used to predict the correlation between base pressure and the 
parameters was determined by a review of the Compendium of Flight Vehicle and 
Base Pressure Techniques (Lockheed Missiles & Space Co., Inc. August 1983). 
From this review, the base pressure technique used was developed. This 
technique is outlined as follows: 

1. Define configuration, dimensions, etc. 

2. Define trajectory and thrust history 

3. Calculate thrust coefficient and thrust loading using base area 

4 . Determine generic base pressure 

5. Determine correction for nozzle axial extension 

6. Determine correction for flare 

7 . Determine correction for multiple nozzle plume effects 

8. Calculate base pressure coefficient. 

The base pressure technique was then applied to each LRB configuration to 
determine the base pressure coefficients. The resulting coefficients were 
compared to actual flight data for validation. The coefficients were then 
used to determine how a change in base configuration would affect aerodynamic 
characteristics and to predict base drag numbers and how base drag would 
affect vehicle performance. Section 7 discusses the effects of the base drag 
further. 


22 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 


LMSC-HEC TR F268592 


The effect of a change in base configuration was analyzed with regard to 
normal force and pitching moment coefficient. The correlation between base 
pressure on the external tank base and pressure distribution on the oribiter’s 
lower wing and fuselage was investigated. The investigation was based on pre- 
vious analysis of solid and gaseous plume simulation test data. Figs. 2-13 
through 2-18 present graphic representations of the analysis. 

From the evaluation of the wind tunnel data the following relationships 
were derived for the LRB study: 


AC = 0.28 AC 

N PB 


4C m . -0.27 4C pB 


ET 


where Ac , Ac are increments relative to the current Shuttle. Ac _ is the 
N M PB^ 

ET base pressure coefficient increment of the LRB configuration relative to the 


current Shuttle. 


4c h 

Ac » 


is the value consistent with plume effects acting in the 


aft region of the vehicle. It was concluded that the mated vehicle normal 
force and pitching moment effects can be predicted if the ET base pressure 
effects can be predicted. 

The LRB plume effects study utilized five different booster configura- 
tions. These configuration are presented in Fig. 2-19. The lack of a definite 
configurations necessitated the estimation of certain parameters. These 
uncertainties are outlined in Table 2-4. 
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Mated Vehicle 

Normal Force Coefficient Increment 
Due to Plume Effects 
(TWT 675) 
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Fig. 2-13 A Normal Force Coefficient (TWT 675) 
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Mated Vehicle 

Pitching Moment Coefficient Increment 
Due to Plume Effects 
(TUT 675) 


□ - X/DSRB - 2 
o - X/DSRB - 3 
^ - X/DSRB - 4 


o 


A 
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O 


Wherei X - distance from ORB base to SSME Exit 
DSRB * SRB base diameter 
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o.oo 


0.05 
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Fig. 2-14 A Pitching Moment Coefficient (TWT 675) 
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Fig. 2-15 A Normal Force Coefficient (IA300) 
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Fig. 2-16 A Pitching Moment Coefficient (IA300) 
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Fig. 2-17 A Normal Force Coefficient (TWT 675 & IA300) 
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Fig. 2 19 Shuttle LRB Configurations Used In Study- 
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Table 2-A ANALYTICAL UNCERTAINTIES 


PARAMETER 

CONFIGURATION 
2 LRB 3 LRB 4 LRB 

5 LRB 

SKIRT DIMENSION 

E 

E 

E 

E 

NOZZLE AXIAL POSITION 

E 

E 

E 

E 

SKIRT FLARE EFFECTS 

EX 

EX 


EX 

ENGINE THRUST 


E 



ENGINE CHAMBER PRESSURE 


E 



NOZZLE AREA RATIO 


E 


E 

PLUME ANGLE VS ALTITUDE 


E 


E 

MULTIPLE ENGINE PLUME EFFECTS 



EX 


BASE PRESSURE ABOVE M ■ 2.0 



EX 


ENGINE SPACING 

E 

E 

E 

E 

ALTITUDE AND DYN PRESSURE VS MACH NO. 

(1) 

(1) 

(1) 

(1) 

THRUST VS MACH NO. 

(2) 

(2) 

(2) 

(2) 


E = ESTIMATED VALUE BASED ON MEASUREMENTS/CALCULATIONS 
Ex * EXTRAPOLATION BEYOND EXISTING DATA BASE WAS PERFORMED. 

(1) = ASSUMED SAME AS CURRENT SHUTTLE PROFILE 

(2) * ASSUMED CONSTANT THRUST 
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The base pressure prediction technique was used to determine the base 
pressure coefficient for each booster configuration. Figures 2-20 and 2-21 
present the Booster base pressure coefficients and external tank base pressure 
coefficients for each configuration at the Mach numbers used in the analysis. 
The base drag increment for the mated vehicle was determined and compared to 
current Shuttle data. Figure 2-22 presents the base drag increments for the 
booster configurations in comparison to the current STS. It was determined 
that the drag for the LRB is greater than that of the current shuttle because 
of the larger base area and the base pressure. 

The normal force and pitching moment effects were calculated and compared 
to the current STS. Figs. 2-23 and 2-24 graphically present the comparison. 
Results showed that, for the LRB's, at M > 1, the pitching moment was sig- 
nificantly increased while normal force was significantly decreased. This was 
attributed to the decrease in ET base pressure due to a larger base area in 
the recirculation base from environment. Furthermore, the greater nozzle area 
ratios of the LRB configurations resulted in lower plume expansion angles and 
decreased base pressure. 

The study recommended that delta base values not be incremented to account 
for plume effects for the following reasons: 


1. Configuration uncertainties and assumptions are a significant factor 

2. Drag increment is not large compared to total vehicle drag and thrust 

3. Normal force increment is not larger compared to total vehicle values 

4. Pitching moment has the most significant impact and could be an 
important factor. 


It is recommended that an update to the plume analysis be conducted when more 
definition becomes available on LRB designs. This updated analysis was 
conducted and is discussed in section 7. 
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I 

I 

I 
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Fig. 2-20 Booster CPB vs Mach 
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Fig. 2-21 ET CPB vs Mach 
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Mated Vehicle Base 
Drag Increment 


CD 

O 



Fig. 2-22 CD Increment vs Mach 
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Fig. 2-23 CN Increment vs Mach 
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2.5 DATA ANALYSIS REVIEW 

2.5.1 Data Analysis Review of Longitudinal Effects 

Test data were analyzed to determine longitudinal effect at negative 

angles of attack. From the test data, it was concluded that the axial force 

coefficient (C ) increases with the diameter and length. The maximum C 

AF AF 

occurs between 170 and 180 ft. See Fig. 2-25 for a plot of C vs booster 

AF 

length and diameter. 

The normal force coefficient (C^) is relatively unaffected by changes 
in length. However, increases in diameter do produce a decrease or more neg- 
ative C^. This decrease is partially due to the larger nose of the vehicle. 

Since the nose of the LRB generated the majority of the normal force, 
increasing the nose area will increase the normal force produced. Whereas 
increasing the length will have very minimal effect on C^. The increased 
plan form area is also a contributing factor. See Fig. 2-26 for a plot of 

C vs booster length and diameter. 

NF 

The pitching moment coefficient decreases with an increase in length due 
to the forward movement of the LRB nose. An increase in diameter has minimal 

effects on C . The LRB length also changes the moment arm. See Fig. 2-27 

M 

for a plot of CMF vs booster length and diameter. 

The aerodynamic center location moves forward with increases in length 
and diameter due to the increased loading and forward movement of the LRB 

nose. See Fig. 2-28 for a plot of X vs booster length and diameter. A 

AC 

summary of vehicle longitudinal effects can be found in Table 2-4. 

2.5.2 Data Analysis Review on Wing Loads 

Analysis of the data with regard to wing loads indicated the wing root coef- 
ficients and elevon hinge moments are relatively unchanged by increase in length. 
However, an increase in diameter increases all wing loads including wing shear, 
root bending, root torsion, and inboard/ outboard elevon hinge moments. 
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TWT0707 ALPHA = -4 



— □ — DILI 
O 01L2 
— A — 01L3 
— V- - 01L4 

— O — 03L1 

D3L2 
— •— 03L3 
A 03L4 


Fig* 2-25 CAF vs Mach (Length and Diameter Effects) 
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TWT0707 ALPHA = -4 
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O D1L2 
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Fig. 2-26 CNF vs Mach (Length and Diameter Effects) 
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Fig. 2-28b X p vs Mach 
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Table 2-5 SUMMARY OF VEHICLE LONGITUDINAL EFFECTS 


• AXIAL FORCE 

- Increases with Diameter (increased frontal area) and Length 

- Maximum C Occurs at LRB Lengths between 170 and 180 ft. 

AF 

0 NORMAL FORCE 

- C^ Relatively Unaffected by Changes in Length 

- Decreases (becomes more negative) with Increases in Diameter 

• Larger Nose is Partially Responsible 

• Increased Plan Form Area a Possible Contributor 

- C Increases with Increase in Diameter 

Na 

- C Relatively unaffected by Change in Length 

• PITCHING MOMENT 

- C^ Decreases with Increase in Length 

• Caused by Forward Movement of LRB Nose 

- Increased Diameter Affects C^ as a Function of LRB Length 

• LRB Length changes Moment Arm 

' c Mec Decrease (slope becomes less negative) with Increases in 
Length and Diameter 

• AERODYNAMIC CENTER LOCATION 

■ X AC moves Forw ard with Increases in Diameter 

• Increased loading on LRB Nose 

-X AC moves Forward with Increases in Length 

• Forward movement of LRB Nose 
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TWT0707 ALPHA = -4 



— □— DILI 
O D1L2 
- -A- - D1L3 
— V — D1L4 
-- 0 -- D4L1 
D4L2 
— ♦ — D4L3 


Fig. 2-29 C gR vs Mach (Length and Diameter Effects) 
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Fig. 2-31 C vs Mach (Length and Diameter Effects) 
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TWT0707 ALPHA * -4 
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Fig. 2-32 CHE I vs Mach (Length and Diameter Effects) 
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TWT0707 ALPHA = -4 
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Fig. 2-33 CHEO vs Mach (Length and Diameter Effects) 
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3. PHASE II LIQUID ROCKET BOOSTER TESTING 

The analysis of Phase I data indicated that increasing the diameter over 
15 ft produced wing load levels that were unacceptable. In order to reduce 
the level of wing loads, various innovative configurations were designed and 
tested. Phase II data were the result of the wind tunnel testing for the new 
configurations. This section delineates the results of the Phase II configura- 
tions and a comparison between Phase I and II data. 

3.1 TEST CONDITIONS AND CONFIGURATIONS 

During Phase II a number of configurations were tested which included 
comparable Phase I designs (Fig. 3-la), Hammerhead designs (Fig. 3-lb), a 
stacked booster design (Fig. 3-lc) and a rotated stack design Fig. 3-ld). 

Test conditions for these configurations can be found in Table 3-1. 

3.2 TEST RESULTS 

The majority of this section presents a comparison of Phase I and Phase 
II results. Therefore, most of the discussion will pertain to config- 
urations SD12L1 and SD15. For a look at the results dealing with other config- 
urations tested see Figs. 3-8a to 3-8i for increments vs Mach number and Figs. 
3- 9a to 3-9e for increments vs LRB diameter. Table 3-3 summarizes the effects 
of LRB diameter on increments for the hammerhead, stacked, and rotated stack 
configurations. 

3.3 PHASE I AND II COMPARISONS 

Comparison between Phase I and Phase II LRB wind tunnel test data revealed 
discrepancies in the incremental data of identical configurations at transonic 
Mach numbers. Through analysis it was determined that the cause of the differ- 
ences was bad baseline data obtained during Phase I testing. 
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Fig. 3-la Phase II Configurations (Phase I Comparison) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



Dx = 12.2 ft 
d 2 = 15.0, 18.0 ft 
L = 159 ft 


Fig. 3-lb Phase II Configuration (MDXXXX Runs) 
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18SO.S23 — j 

X, 21 73.023 


D = 12.2 ft 
L = 159 ft 


Fig. 3-lc Phase II Configuration (STXXXX Runs) 
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D = 15 ft 

0 = 3.6, 7.0, 10.0 deg. 


Fig. 3-ld Phase II Configuration (OSXX Runs) 


55 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 


MATRIX OF DATA FILES USED TO GENERATE PHASE II DATA 


LMSC-HEC TR F268592 


ORIGINAL PAGE IS 
OF POOR QUALITY 


I 

CO 


(0 

H 


M 

U< 

i 

I 

Q 

U 



«o m co 

» •* *— i 


IVM 


l; U « 

£ 2 8 


111 

mu 


co 2 

III 

m 


in ^ 

O CM CM 
*H * 4 H 

> > > 
88 8 


56 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 


This section discusses the analysis performed to determine the validity of 
both sets of data. The conclusions and recommendations resulting from this 
analysis are presented at the end of the section. 

3.3.1 Data Analysis 

Figure 3-2 represents a sketch of the configurations for which data were 
obtained during Phase I tests. Table 3-2 provides the Phase I test matrix for 
wind tunnel tests. In order to distinguish between the Phase I and Phase II 
test data of identical configurations different nomenclatures were used. Sim- 
ilarly, SD12L1 represents the DILI (baseline) configuration. It is important 
to note that all coefficient data have been interpolated to even angles of at- 
tack. No data were generated by interpolating between test Mach numbers. 

Comparisons between Phase I and Phase II incremental data are presented 
in Figs. 3.3a through 3.4f for the D2L2/SD15 configuration. The agreement 
between CNF, CMF, CSR, and CBR data is reasonably good except at the transonic 

Mach numbers of 1.10 and 1.25. The CAF and CTR data are in good agreement 

over the entire Mach number range. 

Shown in Figs. 3.4a through 3.4f are comparisons between Phase I and 
Phase II total data obtained for the D2L2/SD15 configuration for Mach numbers 
ranging from 0.6 to 1.46. The agreement between the two sets of data is quite 

good for all coefficients. Figures 3.5a through 3.5f compare Phase I and 

Phase II total data obtained for the DILI/ SD12L1 configuration. All of the 
data at Mach numbers less than 1.10 agree reasonably well. Additionally, CAF 
and CTR data continue to agree well throughout the entire transonic Mach 
number range. At Mach numbers ranging between 1.10 and 1.25, there is 
considerable difference between Phase I and Phase II results for CNF, CMF, 

CSR, and CBR. 

From the data comparisons made thus far, it is clear that the cause of 
differences seen between the Phase I and Phase II incremental data is invalid 
baseline (D1L1/SD12L1) data. The problem now is to determine which set of 
data is incorrect. 
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Fig. 3-4 Phase I/Phase II Total Data Comparison, 
D2L2/SD15 Configuration (Continued) 
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D2L2/SD15 Configuration (Concluded) 
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Fig. 3 5 Phase I/Phase II Total Data Comparison, 
D1L1/SD12L1 Configuration 
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Fig. 3-5 Phase I/Phase II Total Data Comparison, 
D1L1/SD12L1 Configuration (Continued) 
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Fig. 3-5 Phase I/Phase II Total Data Comparison, 
D1L1/SD12L1 Configuration (Concluded) 
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3.3.2 Angle of Attack Trends . 

From the analyses conducted thus far, it is obvious that the Phase 1 
D1L1/SD12L1 CNF, CMF , CSR, and CBR data are incorrect at Mach numbers between 
1.10 and 1.25 at an angle of an attack of -4 deg. Presented in Figs. 3-6a 
through 3-6h are plots of total coefficient data which show both angle of 
attack trends and diameter trends for different LRB configurations. From 
Figs. 3-6a through 3-6h, it is quite clear that there are significant 
differences between Phase I and Phase II results for the D1L1/SD12L1 
configuration at angles of attack ranging from -4 deg to 0 deg. Although 
positive angle of attack data are not presented in Figs. 3-6a through 3-6h, 
previous analyses have shown that significant differences between Phase I and 
Phase II D1L1/SD12L1 data do exist throughout the entire angle-of-attack range. 

Figures 3-6a through 3-6h also show that the Phase I and Phase II 
D2L2/SD15 data are in good agreement at angles of attack between -4 deg and 0 
deg. Once again, previous analyses have shown that this agreement continues 
at positive angles of attack. 

3.3.3 Results of Data Comparison 

It has been determined that the Phase I D1L1/SD12L1 total CNF, CMF, CSR, 
and CBR data are invalid at Mach numbers between 1.10 and 1.25. Thus, all of 
the Phase I incremental data that were generated at these Mach numbers are 
also invalid. 

Most likely, the cause of the bad Phase I data is potential tunnel 
operating errors. Table 3-2 is a matrix of the Phase I test runs. It has 
been shown that runs TWT033 , TWT0037, and TWT041 contain bad data. If indeed 
the potential tunnel errors resulted in the bad data, then it is quite likely 
that all of the "Dl" data at Mach numbers between 1.10 and 1.25 (runs TWT033 
through TWT044) are also bad due to the consecutive order in which the runs 
were conducted. 
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Figures 3-7a through 3-7f contain comparisons of Phase I "Dl" data and 
Phase II D1L1/SD12L1 data. Figures 3-7a through 3-7f show that some of the 
inconsistencies seen in the Phase I D1L1/SD12L1 data at Mach numbers between 
1.10 and 1.25 are also present in the D1L2, D1L3, and D1L4 data. Thus, the 
validity of the Phase I D1L2, D1L3, and D1L4 total data at Mach numbers 
between 1.10 and 1.25 cannot be confirmed. 

It appears that the extent of the bad Phase I data is limited to runs 
TWT033 through TWT044 (i.e., only the "Dl" configurations). Referring back to 
Figs. 3-6a through 3-6h, the trends produced by the larger diameter configu- 
rations appear to follow the initial trends of the valid (Phase II) D1L1/SD12L1 
data and the D2L2/SD15 data. Unfortunately, there are no Phase II data avail- 
able to confirm the validity of the Phase I "D3" and "D4" data. 

It is recommended that a modified Phase I data base be developed. This 
would be a relatively simple task requiring an additional six test runs. The 
additional runs that would be needed are given in Table 3-3. It is further 
recommended that the increment development method documented in Lockheed- 
Huntsville IDC 88FT44 be used in creating the new Phase I LRB incremental data 
base . 


Table 3-3 ADDITIONAL TEST RUNS SUGGESTED TO REDEVELOP PHASE I DATA BASE 


Configuration M = 1.10 M = 1.25 

D1L2 X X 
D1L3 X X 
D1L4 X X 
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Fig. 3-6 Angle of Attack and Diameter Effects 
on Phase I/Phase II Data (Continued) 
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Fig- 3-9a A CAF vs Diameter 
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Fig. 3-9b A CNF vs Diameter 
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Fig. 3-9c A CN-Wing vs Diameter 
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Fig. 3-9d A CSR vs Diameter 
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4. LRB LATERAL/DIRECTIONAL WIND TUNNEL DATA 

The LRB Lateral/Directional incremental data were taken from Wind Tunnel 
Test TWT0716. The test was conducted from April 1988 to June 1988 to 
investigate length and diameter effects on the lateral/directional aerodynamic 
characteristics of the SSLV. Figure 4-1 presents the scope of the test. The 
configurations used in the test are presented in Figure 4-2. 

The aerodynamic increment coefficients were for the total vehicle and 
were generated using the following equation: 

A DXLY = DXLY - DILI 

where: X = either 2 or 3 (15* or 18* diameter) 

Y = either 1,2,3, or 4 (length variations ranging from 149' to 190*). 

The increments can be used to determine the coefficient increment for any LRB 
configuration. New LRB coefficients equal the current SSLV values plus the 
increments . 


The wing data (bending, shear, torsion) increments are for the right 
wing. The elevon data (inboard and outboard) increments are for the left 
wing. All data are for alpha - 0 and have been uniformly shifted by an 
incremental value so that beta = 0. Sign conventions for the launch vehicle, 
wing, and elevon data are shown in Figure 4-3. 
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Fig. 4-2 Test Configurations (TWT 716) 
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Fig. 4-3a Launch Vehicle Sign Convention 
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Fig. 4-3b Launch Vehicle Sign Convention 
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During analysis of the LRB lateral/directional data, some basic trends 
were observed. The trends are for unshifted increments and are given below. 


1 . CY Trends 


- Generally |CYg| increases in diameter. Thus, for a given B, 

| ACY | increases with increases in diameter. 

- For M < 1.05, LRB length variation has no significant effect on 
CY B . 

- For 1.05 < M < 1.80, |CYg| generally decreases with increases 
in LRB length. 

- For M > 1.80, length trends switch so that |CYg| generally 
increases in LRB length. 

- Typically, LRB length and diameter effects on SSLV CYg values are 
less than 15%. 


2. CYM Trends 


- Generally, |CYMg| increases with increases in diameter. Thus, 
for a given B, |ACYM| increases with increases in diameter. 

- Generally, |CYMg| decreases with increases in LRB length. 

-- Typically, length and diameter effects on SSLV CYMg values are 
less than 20%. 


3. CRM Trends 

- Generally, |CRMg| decreases with increases in diameter. Thus, 
for a given B, |ACRM| increases with increase in diameter. 

- LRB length generally has a small effect on ACRM and thus CRMg. 

- Typically, length and diameter effects on SSLV CRMg values are 
less than 20%. 
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4. CSR Trends - Right Wins 

- For all 8 values. ACSR generally increases as LRB diameter 
increases. 

- For B > 0 and a given diameter, ACSR generally increases as B 
increases for all M. 

- LRB length effects on ACSR are small. 


5. CBR Trends - Right Wing 

- For all B values, ACBR generally increases with increases in LRB 
diameter. 

- For fl > 0 and a given diameter, ACBR typically decreases as B 
becomes more negative. 

- LRB length effects on ACBR typically decreases as B becomes more 
negative. 

- LRB length effects on ACBR values are small. 


6 . CTR Trends - Right Wing 

- For 1.80 < M < 3.48, diameter increases typically cause ACTR to 
decrease for all B. For other Mach numbers, ACTR usually increase 
with increases in LRB diameter. 

- No consistent B trends appear to exist for the Mach numbers tested. 

- LRB length effects on ACTR values are small. 


7. Hinge Moment Trends - Left Wing 

- Inboard - ACHEI 

• ACHEI generally increases with increases in diameter except for 
1.8 < M < 2.5, where it decreases as B increases. 

• ACHEI typically decreases as B increases. 

• LRB length effects on ACHEI are small. 

- Outboard - ACHEO 
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• ACHEO generally increases with increases in diameter except for 
1.8 < M < 2.5, where it increases in diameter. 

• No consistent 8 trends are apparent. 

• LRB length effects on ACHEO are small. 

8. Additional Trends 

• B effects on ACNF, ACMF, ACAF are small for all configurations. 

The data from the wind tunnel test were analyzed with regard to yaw angle 
effects of shuttle wing loads. At Mach = 1.96 the incremental wing loads are 
a strong function of both LRB length and yaw angle. Unlike the longitudinal 
data, the LRB diameter was found to have a small effect on wing loads. 

Analysis of the test data for the various configurations leads to the 
conclusion that the aft skirt on DILI configurations does not significantly 
affect wing loads. It was also determined that the MD15 (hammerhead) con- 
figuration greatly reduces the wing loads over the entire yaw angle range 
tested . The maximum incremental loading on a wing occurs when it is on the 
leeward side. 
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Fig. 4-4 CY vs Length 
B 


100 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 



LENGTH, FT 


Fig. 4-6 CYMg VS LENGTH 
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Fig. 4-7 CYM^ vs Diameter 
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Fig. 4-8 CRH vs Length 
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Fig. 4-9 CRM vs Diameter 
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Table 4-1 DIAMETER EFFECTS OF PHASE II CONFIGURATIONS 


• DIAMETER EFFECTS - AN INCREASE IN DIAMETER INCREASES ACNF, 

MOVES X„„ FORWARD, INCREASES BOTH AC -WING AND Ac 

CP N TR 


• LRB ROTATION EFFECTS (OSXXXX RUNS) - ROTATING THE LRBs DOWN- 
WARD HAS A SMALL EFFECT ON AC VALUES, PLACE X OF THE 

MF CP 

INCREMENTAL LOADING ON THE MOON . HAS A SMALL EFFECT ON BOTH 

AC -WING AND AC 
N TR 


• HAMMERHEAD EFFECTS (MDXXXX RUNS) - MULTI-DIAMETER LRBs CREATE 

LARGE CHANGES IN AC VALUES, HAVE A SMALL EFFECT ON AC VAL- 

NF MF 

UES MOVE X FORWARD, CAUSE A DECREASE IN A C -WING AND AC 

CP N TR 


• STACKED CONFIGURATION EFFECTS (STXXXX RUNS) - STACKING LRBs 

CREATE LARGE CHANGES IN Ac VALUES AND AC , PLACES X„„ 

NF MF CP 

AT THE BASE OF THE ORBITER, CAUSES A SMALL INCREASE IN C -WING 

N 


• HAMMER GAP AND FLOW EXPANSION ARE EFFECTIVE IN REDUCING WING 
LOADS . 
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5 . PROTUBERANCE EFFECTS 


Wind tunnel testing was conducted to investigate the aerodynamic 
effects of protuberances on the SSLV. This section details the results of 
protuberance analysis from wind tunnel testing. Three configurations were 
tested: a baseline SRB, an SRB without protuberances, and an SRB with out 

the TEA Cover. The fairing configurations were also varied. Figures 5-1 
and 5-2 present the four configurations analyzed. A Mach number range of 
0.8 to 4.45 was used for testing. Results from the test can be found in 
Figs. 5-3 to 5-20 which depict protuberance and fairing efforts. 


5 . 1 PROTUBERANCE EFFECTS 


Analysis indicated that SRB protuberances have major effects on wing 
loads. A significant increase in the vehicle normal force increment is ex- 
hibited in the presence of SRB protuberances. The majority of the increase 
can be attributed to orbiter wing shear. Further, the SRB/ET aft attach 
ring/ IEA position and geometry have adverse effects on the orbiter ascent 
wing loads. Wind tunnel testing suggested that ascent wing load reduction 
can be accomplished by removal or modification of the IEA/attach ring. 
Additional efforts, however, would be required to study the impact of the 
IEA relocation. Also a trajectory analysis should be considered to 
determine performance or launch probability increase due to wing load 
reduction. Finally, a complete evaluation of the fairing configurations 
should be performed with regard to aerodynamic enhancement. 
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Fig. 5-1 SRB Protuberance Configurations 
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Fig. 5-2 SRB Fairing Configurations 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA=-4 



— □ — BLSRB 
— o — w/0 IEA 
— a — w/0 PROT 


Fig. 5-3 CNF (Protuberance Effects) 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA. -4 



— a — BLSRB 
— o — W/0 IEA 
— a — w/0 PROT 


Fig. 5-4 CMF (Protuberance Effects) 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA=-4 
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Fig. 5-6 CTOR (Protuberance Effects) 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA=-4 
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Fig. 5-7 CRBM (Protuberance Effects) 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA--4 
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Fig. 5-8 CHMIE (Protuberance Effects) 
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SRB PROTUBERANCE EFFECTS 
JUN 1988 ALPHA— 4 
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Fig. 5-9 CHMOE (Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
JUN 1988 ALPHAa-4 



— o — BLSRB 
— o — RAMP 
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Fig. 5-10 CAF (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
JUN 1988 ALPHA=-4 
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Fig. 5-11 CSR (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig. 5-12 CTOR (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig. 5-13 CRBM (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig. 5-14 CHMIE (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
JUN 1988 ALPHA--4 
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Fig. 5-15 CHMOE (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig- 5-16 CSR (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
JUN 1988 ALPHA--4 
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Fig" 5-17 CTOR (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig. 5-18 CRBM (Fairing/Protuberance Effects) 
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SRB FAIRING EFFECTIVENESS 
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Fig* 5-19 CHMIE (Fairing/Protuberance Effects) 
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Fig. 5-20 CHHOE (Fairing/Protuberance Effects) 
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for an identical SSLV configuration less protuberances. Analysis of the SRB 
protuberance data increments can be summarized as follows: 

• The longitudinal and wing total data show uniform shifts between 
BLSRB and SD12L1-I-AS. Corresponding increment plots show data scatter 
only. 

• Hinge moment protuberance total data do not always give consistent 
changes for changes in alpha. A single number will not always give a 
good representation of the protuberance increment. 

• For M = 1.46, the change in normal force at a = -2 is 0.04. The 
corresponding change in shear is 0.0191 (2 x 0.0191 - 0.0382). Hence 
an increase in normal force is mostly due to an increase in wing 
shear. This shows that the data are consistent. 

• All increment data fall within a reasonable band, as shown in the 
maximum and minimum coefficient table. 

• Some Xcp values are questionable at higher Mach numbers (> 2.74). 

• Wing Xcp values are reasonable, moving aft and slightly inward with 
higher Mach numbers. The exception is M = 4.43, where the wing Sep 
values are located within the body. 

It was then determined that LRB configurations that do not utilize the SRB 
type protuberances need to have the protuberance increments subtracted from 
the aerodynamic data base. To completely exclude the effects of SRB prot- 
uberances from a configuration using Phase I LRB data bases, the coef- 
ficients should be subtracted from the corresponding coefficients of that 
configuration. 
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5.2 SRB INCREMENT DATA BASE SUMMARY 


Data analysis resulted in the generation of a protuberance incremental 
data base. The protuberance incremental data were generated by taking the 
difference between the baseline aerodynamic coefficient data of the SSLV with 
the baseline protuberances included and the aerodynamic coefficients obtained 
for an identical SSLV configuration less protuberances. Analysis of the SRB 
protuberance data increments can be summarized as follows: 


o The longitudinal and wing total data show uniform shifts between 

BLSRB and SD12L1+AS. Corresponding increment plots show data scatter 
only. 

o Hinge moment protuberance total data do not always give consistent 

changes for changes in alpha. A single number will not always give a 
good representation of the protuberance increment. 

o For M = 1.46, the change in normal force at a = -2 is 0.04. The 

corresponding change in shear is 0.0191 (2 x 0.0191 - 0.0382). Hence 
an increase in normal force is mostly due to an increase in wing 
shear. This shows that the data are consistent. 

o All increment data fall within a reasonable band, as shown in the 
maximum and minimum coefficient table. 

o Some Xcp values are questionable at higher Mach numbers (> 2.74). 

o Wing Xcp values are reasonable, moving aft and slightly inward with 

higher Mach numbers. The exception is M = 4.43, where the wing Sep 
values are located within the body. 

It was then determined that LRB configurations that do not utilize the SRB 
type protuberances need to have the protuberance increments subtracted from 
the aerodynamic data base. To completely exclude the effects of SRB prot- 
uberances from a configuration using Phase I LRB data bases, the coef- 
ficients should be subtracted from the corresponding coefficients of that 
configuration. 


134 * 


LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 


6. GAP AND AFT SKIRT EFFECTS 


The effects of the LRB/ET separation gap width and the aft skirt on the 

aerodynamics of the SSLV were analyzed. Data were obtained from wind tunnel 

testing on two configurations at Mach speeds from .9 to 1.5. Figures 6-1 and 

6-2 present the configurations analyzed and the corresponding nomenclatures. 

The gap width was varied on both configurations from 12 in to 33 in. Angles 

of attack ranged from -4 to zero in increments of two. Analysis of the gap 

width was conducted with respect to changes in C , C . C , C , and 

NF AF SR BR 

C values. The objectives of the aft skirt analysis were to determine the 
TR 

protuberances, diameter, and length effects op wing loading. 

6.1 GAP EFFECTS SUMMARY 

The data from the wind tunnel were analyzed to determine the effect of 
gap width on the aerodynamic coefficients. Figures 6-3 to 6-17 show results 
for DILI configurations, Figs. 6-18 to 6-32 for the D2L2 configurations. It 
was determined that increasing the gap width causes an increase in C^ 
values and a decrease in C^, values at negative angles of attack and 
transonic Mach numbers. C^ values tend to increase with increasing gap 
size in the transonic and supersonic image. 

At Mach 1.46, increases in the gap width cause a slight increase in C 

SR 

values for both the DILI (baseline) and D2L2 configuration. However, at Mach 

1.25 a slight increase occurs in the C values for the DILI configuration 

SR 

where as a slight decreases occurs in the D2L2 values. At the same Mach 

numbers, it was observed that the C values for both configurations 

BR 

increase as the gap width increases. The C values and elevon hinge 

TR 

moments decrease slightly on both configurations as gap width increases. 
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Fig. 6-3 C vs Hach (DILI) 
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Fig. 6-4 C vs Mach (DILI) 
NF 
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Fig. 6-6 C vs Mach (DILI) 
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Fig. 6-9 CHEO vs Mach (DILI) 
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Standard Gap vs Standard Gap + 21" 


M-1.46 DILI 



Fig. 6-10 C vs Gap (DILI) 

AF 
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Standard Gap vs Standard Gap + 21” 



Fig. 6-11 C^, vs Gap (DILI) 
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Standard Gap vs Standard Gap + 21" 


M - 1 .46 DILI 



Fig. 6-12 C vs Gap (DILI) 
MF 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-13 C gR vs Gap (DILI) 
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Standard Gap vs Standard Gap + 21” 
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Fig. 6-14 C vs Gap (DILI) 
BR 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-15 C TR vs Gap (DILI) 
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Standard Gap vs Standard Gap + 21 M 
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Fig. 6-16 CHEI vs Gap (DILI) 
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Standard Gap vs Standard Gap + 21' 
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Fig. 6-17 CHEO vs Gap (DILI) 
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Fig. 6-19 C^, vs Mach (D2L2) 
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Fig. 6-23 CHEI vs Mach (D2L2) 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-26 C vs Gap (D2L2) 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-27 C vs Gap (D2L2) 
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Standard Gap vs Standard Gap + 21" 


M - 1 .46 D2L2 



Fig, 6-28 C vs Gap (D2L2) 
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Standard Gap vs Standard Gap + 21" 



Fig. 6-29 C TR vs GAP (D2L2) 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-30 C vs GAP (D2L2) 
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Standard Gap vs Standard Gap + 21" 
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Fig. 6-31 CHEI vs GAP (D2L2) 
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Fig. 6-32 CHEO vs GAP (D2L2) 
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6.2 AFT SKIRT EFFECTS SUMMARY 

Data were also obtained from the wind tunnel test to determine the 
effects of aft skirts on aerodynamic coefficients. The two types of 
configurations tested, DILI and D2L2, can be found in Figs. 6-33 and 6-34. 
The test results can be found in Figs. 6-35 to 6-48 for DILI configurations 
and in Figs. 6-49 to 6-58 for D2L2 configurations. 

The conclusion from the analysis of the aft skirt effect is that the 
addition of the aft skirt had little effect on either total vehicle data or 
wind data. The addition of the skirt also had little effect when analyzing 
diameter and length effects. 
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Fig. 6-33 Aft Skirt Configurations (DILI) 
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Fig. 6-35 C vs Mach (DILI) 
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Fig. 6-37 C vs Mach (DILI) 
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Fig. 6-38 C vs Mach (DILI) 
BR 
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Fig. 6-39 CHE I vs Mach (DILI) 
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Fig. 6-41 C vs a (DILI, M = 1.25) 
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Fig. 6-42 C vs o (DILI, M = 1.25) 
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Fig. 6-43 C TR vs a (DILI, M » 1.25) 
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Fig. 6~A4 C vs a (DILI, M = 1.25) 
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Fig. 6-45 C vs a (DILI, M = 1.47) 
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Fig. 6-49 C vs Mach (D2L2) 
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Fig. 6-54 CHEO vs Mach (D2L2) 
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7. LRB BASE DRAG STUDY 


As part of the effort in selecting a booster in the proposed LRB concept 
for the space shuttle program, Lockheed Missiles and Space Co. conducted a 
base drag study on a number of candidate boosters. At the root of the study 
was the development of a computer code which would calculate base drag based 
on the methods found in the Compendium of Flight Vehicle Base Pressure and 
Base D rag Prediction Techniques (Lockheed Missiles and Space Co.; August 
1983. Once the code was operational It was used to obtain base drag estimates 
on two Martin Marietta and two General Dynamics LRB configurations. These 
results were compared with those found in an earlier base drag study, 
conducted by Lockheed. The final results were then selected, and included in 
the shuttle aerodynamic data base. 

7.1 ORIGINAL BASE DRAG ESTIMATES 

A study hereafter was conducted in late 1987 called the baseline study on 
base drag effects for a number of generic LRB configurations, having 1, 2, 3, 
4, and 5 nozzles (see Fig. 7-1) The study used as its basis the same base 
drag compendium mentioned above. The results were calculated for an STS 
vehicle with SRB’s (see Fig. 7-2) and with an average LRB, Fig. 7-3. The 
trajectory used in these calculations can be found in Figs. 7-4 to 7-6. 

Results for total base drag were calculated for each (see Figs. 7-7, 7-8, and 
Tables 7-1, 7-2). To obtain delta base drag values (Fig. 7.9 and Table 7-3), 
the SRB results were subtracted from the LRB results. These delta base drag 
values were the proposed results to be included in the shuttle aerodynamic 
data base. 
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SHUTTLE BOOSTER CONFIGURATIONS IN STUDY 




Fig . 7-1 Baseline Study LRB Configurations 
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Fig. 7-2 Current SRB 
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Fig. 7-3 Baseline Study LRB 
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Fig. 7-4 Baseline Study Altitude vs Mach 
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Fig. 7-5 Baseline Study Q vs Mach 
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Fig. 7-6 Baseline Study Plume Angle vs Mach 
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MACH 

0.60 

0.74 

0.87 

1.01 

1.14 

1.29 

1.48 

1.71 

1.98 

2.26 

2.82 

3.44 

4.03 


Table 7-1 BASELINE CONCEPT (SRB) 


BASELINE CONCEPT (SRB) 


TOTAL BASE DRAG 


ALT(ft) 

Q( psf ) 

LV(lbs) 

9169.00 

383.80 

85508.00 

13188.39 

503.84 

97314.00 

17772.60 

588.80 

100753.00 

22832.80 

644.33 

106169.00 

27418.19 

683.19 

112520.00 

32386.86 

711.47 

69953.00 

37932.33 

724.01 

30910.00 

44451.32 

699.98 

2197.00 

51310.85 

663.64 

-19574.00 

58188.39 

615.26 

-35652.00 

78577.51 

401.67 

-39647.00 

105414.84 

176.63 

-34905.00 

130953.59 

2.05 

-30393.00 


ORB (lbs) 

20696.00 

23465.90 

24465.40 

28484.40 
32424.20 

21842.33 

14765.33 
8564.36 
3390.80 
-535.40 

-2593.00 

-2050.50 

-1534.25 


2SRB( lbs ) 

19973.00 

20358.00 
18761.10 
17515.40 
16575.20 

5638.74 
-3091 .34 
-8530.64 
-11330.00 
-12541.40 
-10327.48 
-6368.16 
-2600.42 


ET(lbs) 

44839.00 

53491.00 
57527.30 
60169.20 
63520.60 
42472.41 
19236.06 

2163.47 

-11635.05 

-22575.80 

-26727.16 

-26487.22 

-26258.89 
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Table 7-2 BASELINE CONCEPT (LRB) 


BASELINE CONCEPT (LRB) 
TOTAL BASE DRAG 


MACH 

ALT( ft) 

0( psf 5 

LV(lbs) 

ORB ( lbs) 

2 LRB ( lbs ) 

ET( lbs) 

0.60 

9169.00 

383.80 

136013.00 

20696.00 

70478.00 

44839.00 

0.74 

13188.40 

503.85 

154065.30 

23465.90 

77108.40 

53491.00 

0.87 

17772.00 

558.80 

159681.00 

24633.00 

77661.00 

57527.00 

1.01 

22832.80 

644.33 

163921.41 

28484.40 

73288.40 

62160.60 

1.14 

27418.20 

683.19 

167289.41 

32424.20 

66598.80 

68288.00 

1.29 

32386.87 

711.48 

120595.48 

21842.33 

48409.20 

50343.94 

1.48 

37932.34 

724.01 

77027.26 

14765.33 

30752.46 

31557.46 

1.71. 

44451.32 

699.98 

46252.27 

8564.36 

15576.88 

22130.64 

1.98 

51310.85 

663.64 

21057.50 

3390.80 

1627.15 

16028.55 

2.26 

58188.40 

615.26 

787.00 

-535.40 

-10060.40 

11382.80 

2.82 

78577.52 

401.68 

-6680.68 

-2593.00 

-18194.00 

14106.32 

3.44 

105414.84 

176.64 

-2712.06 

-2050.50 

-20240.00 

19578.44 

4.03 

130953.59 

2.05 

1064.53 

-1534.25 

-22187.00 

24785.78 
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Table 7-3 BASELINE CONCEPT (LRB) 


BASELINE CONCEPT (LRB) 
DELTA BASE DRAG 


MACH 

ALT(FT) 

0( psf ) 

LV(lbs) 

ORB ( lbs) 

2RB( lbs ) 

ET(lbs) 

0.60 

9169.00 

383.8 

50505.00 

0.00 

50505.00 

0.00 

0.74 

13188.39 

503.8 

56751.30 

0.00 

56750.40 

0.00 

0.87 

17772.60 

588.8 

58928.00 

0.00 

58899.90 

0.00 

1.01 

22832.80 

644.3 

57752.41 

0.00 

55773.00 

1991.40 

1.14 

27418.19 

683.2 

54769.41 

0.00 

50023.60 

4767.40 

1.29 

32386.86 

711.5 

50642.48 

0.00 

42770.46 

7871.53 

1.48 

37932.33 

724.0 

46117.26 

0.00 

33843.80 

12321.40 

1.71 

44451.32 

700.0 

44055.27 

0.00 

24107.52 

19967.17 

1.98 

51310.85 

663.6 

40631.50 

0.00 

12957.15 

27663.60 

2.26 

58188.39 

615.3 

36439.00 

0.00 

2481.00 

33958.60 

2.82 

78577.51 

401.7 

32966.32 

0.00 

-7866 . 52 

40833.48 

3.44 

105414.84 

176.6 

32192.94 

0.00 

-13871.84 

46065.66 

4.03 

130953.59 

2.0 

31457.53 

0.00 

-19586.58 

51044.67 
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Fig. 7-7 Baseline Study SRB Total Base Drag 
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Fig. 7-8 Baseline Study LRB Total Base Drag 
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Fig. 7-9 Baseline Study LRB Delta Base Drag 
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7.2 BASE DRAG CALCULATION CODE 

Prediction of base drag for a given configuration can be a tedious ef- 
fort, considering all the variables associated with it. A faster means of 
calculation was needed. Using the methods found in the Compendium of Flight 
Veh i cle Base Pressure and Base Draft Prediction Techniques . Lockheed developed 
a FORTRAN computer code. It takes into account all the trajectory and geom- 
etry effects found in the compendium, and obtains base drag predictions. See 
Fig. 7-10 for a flow diagram of the code and Fig. 7-11 for sample results. 


7.3 STS FLIGHTS 2, 3, and 5 

In order to obtain delta base drag values using the new prediction code, 
an STS vehicle with SRB's, on a typical STS trajectory, had to be found to 
submit from. This typical case was found by using an average of STS flights 
2, 3, and 5. The typical trajectory used can be found in Figs. 7-12 to 7-14. 
The base drag results appear in Fig. 7-15 and Table 7-4. 


7.4 MARTIN PUMP FED (LRB1) 

The first LRB candidate configuration used in the base drag code was the 
Martin Pump Fed, shown in Fig. 7-16. The trajectory provided by Martin for 
this case can be found in Figs. 7-17 to 7-19. The base drag results obtained 
for this case can be found in Fig. 7-20 and Table 7-5. After subtracting 
these results with those from STS 2, 3, and 5 the delta base drag values are 
obtained (see Fig. 7-21 and Table 7-6). It is important to note the large 
value for delta base drag found in the Mach 2.0 - > 3.5 range. 
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Fig. 7-10 Base Drag Calculation Code Flowchart 
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MACH 

ALT 

Q 

PRESS 

THRUST TOT 
BOOSTER THRUST 
ORB THRUST 
T/A BOOSTER 
T/A ORB 
T/A ET 
PLUME ANGLE 

UNCORRECTED 
NOZZLE EXT 
NON CYLINDER 
FINS 

TOTAL Pb/Pi 

CPB 

CDB 

DFB (vo/NS) 

DFB (wo/NS) 

NOZ SP Pb/Pi 
NOZ SP Pb 
NOZ SP DRAG 

DFB (w/NS) 

DFB (w/NS) 



2.820000 

(FT) 

76703.00 

(PSF) 

382.0900 

( PSF) 

70. 13728 

(LB) 

5780105. 

(LB) 

2162044. 

(LB) 

1456018. 


(RB 

6550.385 

3356.921 

3415.838 

ORB) 

63.61095 

26.83385 


(RB 

ORB 

ET) 

1.828659 

1 . 178031 

1 . 199655 

(RB 

ORB 

ET) 

- 1 . 95 99998E-02 

1.009191 3E-02 

O.OOOOOOOE+OO 

(RB 

ORB 

ET) 

-0. 1060880 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

(RB 

ORB 

ET) 

0.0000000E+00 

-2. 1762000E-02 

O.OOOOOOOE+OO 

(RB 

ORB 

ET) 

1.702971 

1 . 166361 

1 . 199655 

(RB 

ORB 

ET) 

0. 1262819 

2.9885069E-02 

3.586611 5E-02 

(RB 

ORB 

ET) 

-9.2637083E-03 

-3.346241 5E-0 3 

-7.97 7008 8E-0 3 

(RB ORB 
(TOTAL) 

ET) 

-19042.89 

-30681.18 

-3439.341 

-8198.946 

(RB 

ORB 

ET) 

8.000000 

O.OOOOOOOE+OO 

2 . 6 L I960 

( RB 

ORB 

ET) 

561.0982 

O.OOOOOOOE+OO 

183.1957 

(RB 

ORB 

ET) 

-16015.68 

0.0000000E+00 

-18248. 94 

(RB ORB 
(TOTAL) 

ET) 

-51074.26 
-80961 .48 

-3439.341 

-26447.89 


ig. 7-11 Base Drag Calculation Code Output 
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Fig. 7-13 STS 2, 3, 5 Q vs Mach 
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Fig. 7-14 STS 2, 3, 5 Plume Angle vs Mach 
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Table 7-4 STS 2 , 3 , 5 , (SRB) 


TOTAL BASE DRAG 


MACH 

ALT (FT) 

Q(psf ) 

LV(lbs) 

0.60 

9065.00 

383.79 

82508.93 

0.74 

13447.70 

493.69 

92581.52 

0.87 

17979.00 

576.15 

97112.17 

1.01 

23256.00 

626.44 

107046.77 

1.14 

28248.40 

648.84 

105831.92 

1.29 

33622.93 

657.06 

71643.03 

1.48 

39848.53 

648.33 

30772.77 

1.71 

46490.44 

625.93 

242.14 

1.98 

53950.90 

585.98 

-21445.04 

2.26 

61880.59 

525.74 

-35010.25 

2.82 

83994.19 

323.72 

-36647.77 

3.44 

111853.90 

131.95 

-31002.95 

4.03 

138365.56 

2.05 

-25631.26 


ORB ( lbs) 

2RB( lbs) 

ET( lbs) 

16917.90 

20155.25 

45435.77 

18578.35 

21273.24 

52729.92 

19921.74 

21016.70 

56173.72 

22951.75 

21263.98 

62831.03 

25638.99 

17478.22 

62714.70 

17901.90 

9384.42 

44356.70 

11223.04 

-542.03 

20091.76 

6507.47 

-7259.87 

994.54 

1472.69 

-10187.36 

-12730.37 

-2402.45 

-11259.30 

-21348.49 

-3974.51 

-8090.34 

-24582.91 

-3127.16 

-4140.62 

-23735.16 

-2320.82 

-382.02 

-22928.42 


* AVERAGE VALUES FROM STS 2,3,5 FLIGHTS 
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Fig. 7-15 STS 2, 3, 5 Total Base Drag 
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Fig. 7-16 LRB1 (Martin Pump Fed) 
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Fig. 7-17 LRB1 Altitude vs Mach 
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Fig. 7-18 LRB1 Q vs Mach 
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Fig. 7-19 LRB1 Plume Angle vs Mach 
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MACH 


0.60 

0.74 

0.87 

1.01 

1.14 

1.29 

1.48 

1.71 

1.98 

2.26 

2.82 

3.44 

4.03 


Table 7-5 MARTIN PUMP FED (LRB1) 


TOTAL BASE DRAG 


ALT (FT) 

Q(psf ) 

LV(lbs) 

11439.00 

354.1 

90121.27 

15188.00 

470.9 

97449.14 

18884.00 

565.0 

100679.61 

23056.00 

635.9 

126811.89 

27677.00 

672.8 

116658.11 

32728.00 

692.0 

91287.60 

38282.00 

701.0 

50612.79 

44439.00 

697.6 

36313.83 

51272.00 

666.2 

1891.44 

58804.00 

588.9 

-38116.73 

75818.00 

399.1 

-85398.99 

95308.00 

241.4 

-64089.47 

117102.00 

128.1 

-33940.29 


ORB( lbs) 

2RB(lbs) 

ET(lbs) 

22647.26 

25439.72 

42034.27 

24769.23 

27356.17 

45323.73 

26067.21 

28357.47 

46254.92 

33016.43 

34399.14 

59396.33 

30605.34 

30472.93 

55579.84 

25161.46 

20321.63 

45804.51 

16100.30 

5776.70 

28735.79 

13497.04 

-235.69 

23052.48 

8644.59 

-18148.03 

11394.88 

2935.08 

-40194.39 

-857.42 

-3402.95 

-65108.59 

-16887.46 

-3830.44 

-38310.08 

-21948.95 

-2794.35 

-17728.64 

-13417.30 
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Table 


MACH 

ALT (FT) 

0( psf 

0.60 

9065.00 

354.1 

0.74 

13447.70 

470.9 

0.87 

17979.00 

565.0 

1.01 

23256.00 

635.9 

1.14 

28248.40 

672.8 

1.29 

33622.93 

692.0 

1.48 

39848.53 

701.0 

1.71 

46490.44 

697.6 

1.98 

53950.90 

666.2 

2.26 

61880.59 

588.9 

2.82 

83994.19 

399.1 

3.44 

111853.90 

241.4 

4.03 

138365.56 

128.1 


7-6 MARTIN PUMP FED 


DELTA BASE DRAG 
) LV(lbs) ORB (lbs) 


7612.34 5729.36 

4867.62 6190.88 

3567.44 6145.47 

19765.12 10064.68 

10826.19 4966.35 

19644.57 7259.56 

19840.02 4877.26 

36071.69 6989.57 

23336.48 7171.90 

-3106.48 5337.53 

-48751.22 571.56 

-33086.52 -703.28 

-8309.03 -473.53 


(LRB1) 


2RB ( lbs ) 

ET(lbs) 

5284.47 

-3401.50 

6082.93 

-7406.19 

7340.77 

-9918.80 

13135.16 

-3434.70 

12994.71 

-7134.86 

10937.21 

1447.81 

6318.73 

8644.03 

7024.18 

22057.94 

-7960.67 

24125.25 

-28935.09 

20491.07 

-57018.25 

7695.45 

-34169.46 

1786.21 

-17346.62 

9511.12 


DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 
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Fig. 7-20 LRB1 Total Base Drag 
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* DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 


Fig. 7-21 LRB1 Delta Base Drag 
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7.5 MARTIN PRESSURE FED (LRB2) 

The next LRB candidate configuration used in the base drag code was the 
Martin Pressure Fed, shown in Fig. 7-22. The trajectory provided by Martin 
for this case can be found in Figs. 7-23 to 7-25. The base drag results 
obtained for this case can be found in Fig. 7-26 and Table 7-7. After 
subtracting these results with those from STS 2, 3, and 5 the delta base drag 
values are obtained (see Fig. 7-27 and Table 7-8). It is important to note 
the large value for delta base drag found in the Mach 2.0 - > 3.0 range. 

7.6 GENERAL DYNAMICS 02/H2 PUMP FED (LRB3) 

The next LRB candidate configuration used in the base drag code was the 
General Dynamics 02/H2 Pump Fed, shown in Fig. 7-28. The trajectory provided 
by General Dynamics for this case can be found in Figs. 7-29 to 7-31. The 
base drag results obtained for this case can be found in Fig. 7-32 and Table 
7-9. After subtracting these results with those from STS 2, 3, and 5 the 
delta base drag values are obtained (see Fig. 7-33 and Table 7-10). It is 
important to note the large value for delta base drag found in the Mach 1.5 - 
> 3.5 range. 

7.7 GENERAL DYNAMICS 02/RP1 PUMP FED (LRB4) 

The next LRB candidate configuration used in the base drag code was the 
General Dynamics Pump Fed shown in Fig. 7-34. The trajectory provided by 
Martin for this case can be found in Figs. 7-35 to 7-37. The base drag 
results obtained for this case can be found in Fig. 7-38 and Table 7-11. 

After subtracting these results with those from STS 2, 3, and 5 the delta base 
drag values are obtained (see Fig. 7-39 and Table 7-12. It is important to 
note the large value for delta base drag found in the Mach 2.0 - > 3.5 range. 
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Fig. 7- 22 LRB2 (Martin Pressure Fed) 
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MARTIN PRESSURE FED (LRB2) 



Fig . 7-23 LRB2 Altitude vs Mach 
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Fig. 7-24 LRB2 Q vs Mach 
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Fig. 7-25 LRB2 Plume Angle vs Mach 
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Table 7~7 MARTIN PRESSURE FED (LRB2) 


TOTAL BASE DRAG 


MACH 

ALT (FT) 

Q( psf ) 

LV(lbs) 

ORB ( lbs) 

2RB( lbs ) 

ET(lbs) 

0.60 

8224.00 

398.9 

111153.84 

25512.55 

38288.93 

47352.36 

0.73 

12404.00 

503.9 

115526.58 

26845.83 

39516.39 

49164.35 

0.94 

20182.00 

625.4 

127219.05 

30258.93 

43060. 29 

53899.82 

1.07 

24713.00 

671.0 

138130.73 

33027.53 

45859.09 

59244.11 

1.21 

29714.00 

698.0 

118996.87 

29539.59 

36713.43 

52743.85 

1.38 

35229.00 

706.1 

79519.70 

20898.54 

21128.10 

37493.05 

1.57 

41319.00 

689.0 

42197.63 

14026.65 

2531.22 

25639.75 

1.81 

48042.00 

652.5 

16510.84 

10483.26 

-13790.10 

19817.68 

2.05 

55422.00 

576.2 

-9962.88 

5000.08 

-24159.76 

9196.80 

2.28 

63460.00 

473.6 

-49951.26 

276.37 

-49719.78 

-507.85 

2.75 

81442.00 

292.6 

-69105.40 

-3948.87 

-55328.49 

-9828.04 

3.30 

101970.00 

165.5 

-45803.08 

-5598.42 

-29895.52 

-10309.15 

3.86 

125165. '00 

84.4 

-27067.54 

-2439.16 

-15062.38 

-9566.00 
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Table 7-8 MARTIN PRESSURE FED (LRB2) 


DELTA BASE DRAG 


MACH ALT (FT) Q(psf) LV(lbs) 


0.60 9065.00 398.9 28644.91 
0.73 13134.65 503.9 23664.53 
0.94 20460.00 625.4 26567.23 
1.07 25687.20 671.0 25395.14 
1.21 30776.00 698.0 28257.98 
1.38 36752.53 706.1 31916.80 
1.57 42516.68 689.0 24036.60 
1.81 49317.55 652.5 28352.84 
2.05 55858.75 576.2 15436.35 
2.28 62525.13 473.6 -14561.76 
2.75 80848.75 292.6 -31820.31 
3.30 105563.00 165.5 -13525.49 
3.86 130726.59 84.4 111.50 


ORB (lbs) 

2RB( lbs ) 

ET (lbs) 

8594.65 

18133.68 

1916.59 

8386.08 

18323.01 

-3044.56 

8827.58 

23086.47 

-5346.82 

8093.05 

24272.83 

-6970.72 

7358.78 

23106.31 

-2207.09 

7834.00 

16607.88 

7474.92 

3979.94 

6961.17 

13095.48 

6404.33 

-4649.43 

26597.95 

4600.55 

-13541.41 

24377.21 

2846.27 

-38554.68 

21146.64 

121.31 

-46792.22 

14850.58 

-2279.92 

-24863.03 

13617.43 

113.99 

-13597.38 

13594.87 


* DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 
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Fig. 7-26 LRB2 Total Base Drag 
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* DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3.5 


Fig. 7-27 


LRB2 Delta Base Drag 
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Fig. 7-28 LRB3 (General Dynamics 02H2 Pump Fed) 
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Fig. 7-29 LRB3 Altitude vs Mach 


234 


LOCKHEED- HUNTSVILLE ENGINEERING CENTER 




(PSF) 


LMSC-HEC TR F268592 


GENERAL DYNAMICS 02H2 PUMP FED (LRB3) 



0 1 2 3 4 5 

MACH NUMBER 


Fig. 7-30 LRB3 Q vs Mach 
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Fig. 7-31 LRB3 Plume Angle vs Mach 
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MACH 


0.60 

0.74 

0.87 

1.01 

1.14 

1.29 

1.48 

1.71 

1.98 

2.26 

2.82 

3.44 

4.03 


Table 7-9 GENERAL DYNAMICS 02H2 PUMP (LRB3) 


TOTAL BASE DRAG 


ALT (FT) 

0( psf ) 

LV(lbs) 

8967.00 

382.1 

99706.46 

12248.00 

511.6 

108641.82 

15627.00 

618.8 

113187.52 

21627.00 

648.6 

135482.20 

25321.00 

651.0 

115992.51 

33079.00 

637.1 

72657.88 

38191.00 

659.0 

10533.36 

44691.00 

641.8 

-22206 . 24 

52458.00 

594.8 

-58799.78 

60462.00 

528.7 

-85485.80 

76703.00 

382.1 

-80961.48 

94896.00 

245.1 

-64203.41 

113267.00 

146.6 

-35320.96 


ORB ( lbs) 

2RB( lbs) 

ET( lbs) 

24435.51 

29917.63 

45353.33 

26912.14 

32484.75 

49244.93 

28548.44 

33981.34 

50657.73 

33748.65 

39947.76 

61785.79 

29424.86 

31099.15 

55468.50 

21992.71 

7210.25 

43454.92 

13895.38 

-26411.73 

23049.71 

10495.01 

-49562.45 

16861.20 

6412.80 

-72319.77 

7107.20 

2142.89 

-84820.01 

-2808,68 

-3439.34 

-51074.26 

-26447.89 

-3781.56 

-28020.01 

-32401.84 

-2895.41 

-15719.96 

-16705.59 
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MACH 


0.60 

0.74 

0.87 

1.01 

1.14 

1.29 

1.48 

1.71 

1.98 

2.26 

2.82 

3.44 

4.03 


Table 7-10 GENERAL DYNAMICS 02H2 PUMP (LRB3 ) 


DELIA BASE DRAG 

ALT(FT) Q(psf) LV(lbs) ORB(lbs) 2RB(lbs) ET(lbs) 


9065.00 

382.1 

17197.53 

13447.70 

511.6 

16060.30 

17979.00 

618.8 

16075.35 

23256.00 

648.6 

28435.43 

28248.40 

651.0 

10160.59 

33622.93 

637.1 

1014.85 

39848.53 

659.0 

-20239.41 

46490.44 

641.8 

-22448.38 

53950.90 

594.8 

-37354.74 

61880.59 

528.7 

-50475.55 

83994.19 

382.1 

-44313.71 

111853.90 

245.1 

-33200.46 

138365.56 

146.6 

-9689.70 


7517.61 

9762.38 

-82.44 

8333.79 

11211.51 

-3484.99 

8626.70 

12964.64 

-5515.99 

10796.90 

18683.78 

-1045.24 

3785.87 

13620.93 

-7246.20 

4090.81 

-2174.17 

-901.78 

2672.34 

-25869.70 

2957.95 

3987.54 

-42302.58 

15866.66 

4940.11 

-62132.41 

19837.57 

4545.34 

-73560.71 

18539.81 

535.17 

-42983.92 

-1864.98 

-654.40 

-23879.39 

-8666.68 

-574.59 

-15337.94 

6222.83 


DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 
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Fig. 7-32 LRB3 Total Base Drag 
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* DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 


Fig. 7-33 LRB3 Delta Base Drag 
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Fig. 7-34 LRB4 (General Dynamics 02RP1 Pump Fed) 
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Fig. 7-35 LRB4 Altitude vs Mach 
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Fig. 7-36 LRB4 Q vs Mach 


243 


LOCKHEED- HUNTSVILLE ENGINEERING CENTER 




PLUME ANGLE 


LMSC-HEC TR F268592 


GENERAL DYNAMICS 02RP1 PUMP FED (LRB4) 



Fig. 7-37 LRB4 Plume Angle vs Mach 
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Table 7-11 GENERAL DYNAMICS 02RP1 PUMP FED (LRB4) 


TOTAL BASE DRAG 


MACH 

ALT (FT) 

0( psf ) 

LV(lbs) 

ORB ( lbs) 

2RB ( lbs ) 

ET(lbs) 

0.60 

8825.00 

384.2 

111958.64 

24570.46 

41784.38 

45603.80 

0.74 

12268.00 

511.3 

120945.38 

26895.31 

44835.93 

49214.13 

0.87 

15612.00 

617.2 

125708.64 

28476.93 

46700.87 

50530.83 

1.01 

19749.00 

703.2 

163307.53 

36849.52 

59631.35 

66826.65 

1.14 

25233.00 

707.0 

143871.08 

32579.33 

51236.24 

60055.51 

1.29 

30997.00 

700.3 

113054.79 

25700.65 

37776.69 

49577.44 

1.48 

37812.00 

670.9 

48943.43 

14535.58 

3492.97 

30914.87 

1.71 

46603.00 

588.5 

24445.61 

7992.53 

-3005.13 

19458.21 

1.98 

55905.00 

505.2 

-18465.70 

3344.40 

-33081.82 

11271.73 

2.26 

65080.00 

425.0 

-54432.28 

-445.97 

-57637.99 

3651.69 

2.82 

79785.00 

329.8 

-73499.05 

-3327.81 

-55249.38 

-14921.86 

3.44 

96629.00 

226.9 

-47778.41 

-2166.37 

-31608.99 

-14003.06 

4.03 

116507.00 

129.2 

-31993.45 

-1015.14 

-18291.99 

-12686.32 
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MACH 


0.60 

0.74 

0.87 

1.01 

1.14 

1.29 

1.48 

1.71 

1.98 

2.26 

2.82 

3.44 

4.03 


Table 7-12 GENERAL DYNAMICS 02RP1 PUMP FED (LRB4) 


DELTA BASE DRAG 


ALT (FT) 

Q(psf ) 

9065.00 

384.2 

13447.70 

511.3 

17979.00 

617.2 

23256.00 

703.2 

28248.40 

707.0 

33622.93 

700.3 

39848.53 

670.9 

46490.44 

588.5 

53950.90 

505.2 

61880.59 

425.0 

83994.19 

329.8 

111853.90 

226.9 

138365.56 

129.2 


LV(lbs) 

ORB(lbs) 

29449.71 

7652.56 

28363.86 

8316.96 

28596.47 

8555.19 

56260.76 

13897.77 

38039.16 

6940.34 

41411.76 

7798.75 

18170.66 

3312.54 

24203.47 

1485.06 

2979.34 

1871.71 

-19422.03 

1956.48 

-36851.28 

646.70 

-16775.46 

960.79 

-6362.19 

1305.68 


2RB( lbs ) 

ET(lbs) 

21629.13 

168.03 

23562.69 

-3515.79 

25684.17 

-5642.89 

38367.37 

3995.62 

33758.02 

-2659.19 

28392.27 

5220.74 

4035.00 

10823.11 

4254.74 

18463.67 

-22894.46 

24002.10 

-46378.69 

25000.18 

-47159.04 

9661.05 

-27468.37 

9732.10 

-17909.97 

10242.10 


DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 
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Fig. 7-38 LRB4 Total Base Drag 
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* DELTA DRAG VALUES ARE BASED ON DELTA FROM STS FLIGHTS 2,3,5 


Fig. 7-39 LRB4 Delta Base Drag 
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7.8 CONCEPT COMP ARISON /CONCLUSION 

The next step was to choose which of the base drag values to use in the 
shuttle aerod^amic data base. The baseline study provided average values for 
a set of LRB configurations, and the LRB base drag study a different set of 
values for each specific configuration. If the base drag values have little 
variance from one configuration to another and are a small percentage of the 
total Launch vehicle base drag, the baseline study results will be sufficient. 
If the opposite case is true the values from the LRB base drag study should 
be used. Shown in Figs. 7-40 to 7-43 are a comparison of total base drag from 
both studies. In Fig. 7-42 the same 100,000 lb variation occurs in booster 
base drag. In Fig. 7-43 up to 40,000 lb variation is found when comparing ET 
base drag. 

When comparing the results from both studies a large variation can be 
seen from one configuration to another and the difference is a significant 
part of the total base drag. It is for this reason that Lockheed recommends 
using the results from the new LRB base drag study. 
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Fig. 7-40 Launch Vehicle Total Base Drag 
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Fig. 7 41 Orbiter Base Total Base Drag 
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Fig. 7-42 Booster Total Base Drag 


252 

LOCKHEED-HUNTSVILLE ENGINEERING CENTER 



TOTAL BASE DRAG (LBS) 


LMSC-HEC TR F268592 


EXTERNAL TANKS 



MACH NUMBER 


-GI- 


ST S235 

9LRB 



MLRB1 

MLRB2 

GDLR83 

G0LRB4 


Fig. 7-43 External Tank Total Base Drag 
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8 . SUMMARY 


Phase I and Phase II of wind tunnel tests comprise a number of LRB config- 
urations varying in length and diameter, with and without protuberances and 
aft skirts and varying ET-LRB gap width. 


Conclusions drawn from these tests, with regard to varying length/dia 
meter, included the following Longitudinal Trends: 

• CAF increased with both length and diameter. 

• CNF was relatively unaffected by changes in length, but a more negative 
CNF was produced by an increase in diameter. 

e CMF was relatively unaffected by changes in diameter, but CMF decreases 
with an increase in length. 

e The aerodynamic center moves forward with an increase in both length 
and diameter. 

e Wing loading coefficients CSR, CRB, CTR, and elevon coefficients CHE1, 
and CHEO were relatively unchanged by increase in length but were all 
increased by increase in diameter. 

Conclusions drawn from the tests, with regard to varying length/diameter, 
included the following lateral/directional trends: 


• | CYB | increases with increase in diameter. Length has little 
effect on |ACYB| below M = 1.05, |CYB| decreases with 

increase in length for 1.05 < M < 1.80, and increases with increase n 
length for M > 1.80. 

• |CYMB| increases with increase in diameter, and |ACYMF| 
decreases with increase in length. 

• | CRMB | decreases with increase in diameter, and is relatively unaf- 
fected by increase in length. 

• | ACSR | generally increases with increases in diameter (for all 
B), and has little length effect. 
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• | ACTR| for 1.80 < M < 3.48, an increase in diameter produces a 
decrease in |ACTR| (for all B) ; for all other M an increase in 
diameter shows an increase in CTR. Length effects on CTR are small. 

• ACHEI general 1 / increases with increase in diameter except for 1.8 < 

M < 2.5 where ACHEI decreases. Length has little effect. 

t ACHEO generally increases with increase in diameter except for 1.8 < 

M < 2.5 where ACHEO decreases. Length has little effect. 

Analysis of SRB protuberances showed that they had a major effect on both 
wing loads and SSLV normal force coefficient. The test first recommended 
modification of certain protuberances. It was later decided that since 
proposed LRB configurations did not have these protuberances, the coefficients 
should be removed for the LRB data base. 

Conclusions drawn from the tests, with regard to varying gap width, 
included the following trends: 


• CAF increases as gap width increases for transonic and supersonic Mach 
numbers . 

• CNF increases as gap width increases in the transonic range (for 
negative alpha) . 

• CMF decreases as gap width increases in the transonic range (for 
negative alpha). 

• Effects of gap width on CSR vary with configuration. 

• CBR increases as gap width increases . 

• CTR, CHEI, CHEO decreases as gap width increases. 


Conclusions drawn from these tests, with regard to aft skirt variation, 
showed that, with the exception of a slight increase in CAF, the aft skirt had 
little effect on either wing or total vehicle data. 

Finally, on the subject of base drag/plume effects it was concluded that 
the effects were significant. It was recommended that each specific LRB 
configuration's base drag be calculated before adding to the LRB data base. 
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BASE DRAG CALCULATION CODE USERS GUIDE 


A-l INTRODUCTION 

The Base Drag calculation code, BASE4, is based on the methods found in 
the Compendium of Flight Vehicle Base Pressure and Base Drag Prediction 
Techniques, (Lockheed Missiles and Space Co., 1983). Sections 2 and 3 of this 
document give a detailed background of base drag theory, which the code 
follows. The remainder of this appendix will deal only with issues involved 
in running the BASE4 Code, and assumes the user has some knowledge of the 
parameters involved. If any questions arises the reader should refer to the 
above document for more detail. 

A- 2 CODE ALGORITHM 

BASE4 is presently set up to handle a SSLV, with the Orbiter, external 
tank and a user selected pair of boosters. The boosters can vary in size, 
configuration, and thrust profile. The user supplies the booster geometry, 
thrust profile, and trajectory. Among the factors the code takes into account 
during calculation are corrections for nozzle extension, non-cylindrical 
shape, addition of fins, and nozzle spacing. A flow chart for BASE4 can be 
found in Fig. A-l. A program listing can be found at the end of this section. 

A- 3 PARAMETERS 

A list of parameters used in the BASE4 code can be found in Table A.l. 

All units in the code assumed to be in the English system except where noted 
otherwise. Those parameters in Table A-l noted with an asterisk are 
inputs supplied by the user, those with a "**" are tabular inputs supplied by 
the user. Figures A-2 and A- 3 depict these input parameters and their usage. 


A-l 


LOCKHEED- HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 


A-4 PROGRAM INPUT /OUTPUT 

All I/O in BASE4 is handled via files. For input the user creates and 

specifies the input file, which the code prompts for, and the code generates 

five ouput files. The following describes these I/O files: 

Filename Format 

xxxxx.INP User supplied input file. See Fig. 1-4 for a sample input 

case. 

BASE4 . OUT Tabular listing of Mach, Altitude, Q, Total Base Drag, 

Orb iter Base Drag, Booster Base Drag, and External Tank Base 
Drag at each of the specified trajectory points. See Fig. 

A-5 for a sample case. 

BASE4 . TRACE Listing of the calculation of base drag at each of the 

specified trajectory points. See Fig. A-5 for a sample case. 

BASE 4 . TRAJ Tabular listing of Mach, Altitude, and Dynamic Pressure. See 

Fig. A-7 for a sample case. 

BASE4. THRUST Tabular listing of Mach, Alttitude, Orbiter Thrust, Single 

Booster Thrust, and Total Thrust, at each of the specified 
trajectory points. See Fig. A-8 for a sample Case. 

BASE4. PLUME Tabular listing of Mach, Altitude, Orbiter Plume Angle and 

Booster Plume Angle at each of the specified trajectory 
points. See Fig. A-9 for a sample case. 
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TABLE A . 1 LIST OF PARAMETERS (cont. 


LMSC HEC TR F268592 
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Table A-l LIST OF PARAMETERS (Continued) 



LMSC-HEC TR F268592 



Fig. A-l BASE4 Flowchart 
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Fig. A- 2 Geometric Inputs 
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INPUTS FOR MOTOR CHARACTERISTICS 


EELRB - Nozzle Expansion Ratio 
PCLRB - Motor Chamber Pressure 
KOLRB - Motor Constant 


(All three are inputs which determine the Booster's plume angle 
effects at each given altitude in the trajectory) 

Dj = KO * (alt) * (Pc) 0 * 8 * (EE) -0 * 5 

The inputs EEET, PCET and KOET, are the similar inputs for the 
plume effects on the External Tank. These Values in general use 
the same EE, Pc and KO as the Booster. 


AFLRB - Base Area Multiplying Factor 

For l Nozzle - AFLRB = 1.0 

For 2 Nozzles - AFLRB = 0.333 

For 3 Nozzles - AFLRB = 0.333 

FOR 4 Nozzles - AFLRB = 0.333 

(This quantity accounts for the CP distribution that occurs in 
the inner area between multi-nozzle set-ups) 


TABULAR TRAJECTORY INPUTS 

Trajectory inputs should be in the following form: 

Mach(l) Altitude ( 1 ) Q ( 1 ) THTOT(l) THLRB( 1 ) 

• * • • « 

• • • • « 

• • • • • 

Mach(n) Altitude(n) Q(n) THTOT(n) THLRB(n) 

THTOT is the total mated vehicle thrust 
THLRB is the thrust of one booster 

Fig. A-3 Motor Characteristics and Trajectory Inputs 


A- 7 


LOCKHEED- HUNTSVILLE ENGINEERING CENTER 



LMSC-HEC TR F268592 


ORIGINAL PAGE 13 
OF POOR QUALITY 


* INPUT FILE FOR BASE4 

* 

* VEHICLE - SSLV 

* 

** GEOMETRIC PARAMETERS 
* 


1 

17.3 

12.2 

15.3 
4.9 
0.0 

0.0001713 
7.16 
105840 . 
0.000171 3 
7.16 
105840. 
1.0 

* 

** 


INLRB 

DBLRB 

DELRB 

DCYLRB 

XJLRB 

XSLRB 

KOLRB 

EELRB 

PCLRB 

KOET 

EEET 

PCET 

AFLRB 

PARAMETERS 


TRAJECTORY 


** MACH, 
* 

ALTITUDE, Q, 

TOTAL-THRUST 

) . 6 0 

9065.00000 

383. 79999 

j . 7 4 

13447. 70020 

493.69998 

0.8 7 

17979.00000 

576. 15002 

l .01 

23256.00195 

626.44000 

L . 14 

28248.40039 

648.84003 

l . 29 

33622.93359 

657.06665 

1.48 

39848. 53516 

648. 33997 

1.71 

46490.44141 

625.93597 

l .98 

53950.90234 

585.98499 

2.2 6 

61880.59766 

525.74005 

2.82 

83994. 19531 

323.72403 

1 . 4 4 

1 1 1853.90625 

1 3 1 .95799 

i . J 3 

i. 3 8 365.56250 

' . 0 5000 


BOOSTER-THRUST 

6437333.00000 

6039033.00000 

5786666.00000 

5569533.00000 

5605800.50000 

5926244.50000 

6219244. 50000 

6416799.50000 

6494199.50000 

6473133.00000 

5825667.00000 

5030000.00000 

4272832. 50000 


2607000 . 00000 

2474000.00000 

2369633.00000 
956333. 18750 

2283133. 50000 

2369200.00000 

2461622. 50000 
2522293.25000 

2554483.00000 
2538999.75000 

2211053.50000 

1810740.00000 
1429796. 3 7 500 


Fig. A-4 Sample Input File 
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LMSC-HEC TR F268592 


MACH 

0. AO 
0,7 4 
0.87 
l • • 1 

' . ! 4 
1.29 

I . 48 

l • 7 1 
I . 98 

2. 26 
2.82 
9.4 


ORIGINAL PAGE IS 
OF POOR QUALITY 


XXXXXXXX XXXXXXXX (XXX) 
TOTAL BASE DRAG 


ALT(FT) 

Q(psF) 

LV(lbs) 

9065.00 

383.8 

88103. 11 

13447.70 

493.7 

92256.7 1 

17979.00 

576.2 

92559.63 

23256.00 

626.4 

106989.99 

28248.40 

648.8 

95184. 19 

33622.93 

657 .1 

66081.19 

39848. 54 

648.3 

24157.26 

46490.44 

625.9 

4889.39 

53950.90 

586.0 

-18516.41 

61880.60 

525.7 

-35989.62 

83994.20 

323.7 

-39757.47 

l 1 1853.91 

132.0 

-26571.28 

: 38 36 5 . 56 

2 . 0 

-3863.47 


ORB(lbs) 

2RB( lbs ) 

ET(lbs) 

24546.79 

17996.44 

45559.88 

25968.50 

18769.99 

47518.21 

26581 . 64 

18810.24 

47167.74 

26947.87 

22794.83 

57247.29 

31010.86 

13959.38 

50213.95 

25314.08 

4729.38 

36037.72 

14983.85 

-5636.19 

14809.60 

10445.44 

-11735.59 

6179.54 

6899.03 

-16862. 33 

-8553. 1 1 

2291 . 34 

-16975.81 

-21305. 15 

-3714.92 

-1 3536. 1 7 

-22506. 38 

-3002.27 

-6822.41 

-16746.60 

- 5 8 3 5 

-144.85 

-3660.27 


Fig. A- 5 Sample BASE4.0UT File 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


LMSC-HEC TR F268592 


\ ft l % B 

235.0618 

ABORB 

433.736 1 

ABET 

598. 2849 

Ah LRB 

116.8987 

AEORB 

132.5359 

AEFET 

1502.145 

AEFLRB 

118. 1631 

AEFORB 

301.2002 

AEET 

350.6960 

AO Y LRB 

183.8539 

L 1 ORB 

17.50000 

DEET 

21.13102 

ARLRB 

1 . 278525 

L20RB 

17.50000 

ANET 

299. 1425 

1, 1 LRB 

12.20000 

L30RB 

15. 15544 

XSDET 

1.975768 

L 2 L RB 

12.20000 

ANORB 

7 . 420349 IE-02 



AN LRB 

31.94133 

XSDORB 

1.333333 



XSDLRB 

0.0000000E+00 






**** *************************************** 
A****************************************** 


MACH 



0.6000000 



ALT 

(FT) 

9065 

. 000 



0 

(PSF) 

383. 

8000 



PRESS 

(PSF) 

1535 

. 280 



THRUST TOT 

(LB) 

6437333. 



BOOSTER THRUST 

' (LB) 

2607000. 



'8 THRUST 

(LB) 

1223333. 



C f 8 00STER 


28. 

89708 




CTORB 


7.348761 




c T (•: T 


1 1 . 

16578 




POJ'ME ANGLE 

(RB 

ORB) 


6.072794 

3.171308 


UNCORRECTED 

(RB 

ORB 

ET) 

0.9500000 

0.9500000 

0.9500000 

M 0 • L L E EXT 

(RB 

ORB 

ET) 

0.0000000E+00 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

NON CYLINDER 

(RB 

ORB 

ET) 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

FINS 

(RB 

ORB 

ET) 

0.0000000E+00 

-3.5 1 0000 2E-0 3 

O.OOOOOOOE+OO 

TOTAL Pb/Pi 

(RB 

ORB 

ET) 

0.9500000 

0.9464900 

0.9500000 

CPB 

(RB 

ORB 

ET) 

-0.1984127 

-0.2123413 

-0. 1984127 

CDB 

(RB 

ORB 

ET) 

8. 7 156398E-03 

2. 3775930E-02 

4.41291 26E-0 2 

OFB (wo/NS) 

(RB 

ORB 

ET) 

17996.44 

24546.79 

45559.88 

OFb (wo/NS) 

(TOTAL) 


88103. 1 1 



NOZ SP Pb/Pi 

(RB 

ORB 

ET) 

0.0000000E+00 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

NOZ SP Pb 

(RB 

ORB 

ET) 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

NO?: SP DRAG 

(RB 

ORB 

ET) 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

!)FB (w/NS) 

(RB 

ORB 

ET) 

17996.44 

24546.79 

45559.88 

Dl-B (w/NS) 

(TOTAL) 


88103.11 




******************************************* 


MACH (n) 

l ' V i ! . ; n ) 
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^ ORIGINAL PAGE IS 

j 0F POOR QUALITY 

LMSC-HEC TR F268592 

MACH 

ALTITUDE 

(FT) 

Q 

( PSF ) 

0.6000000 

9065.000 

383.8000 

g 0.7400000 

13447.70 

493.7000 

■ 0.8700000 

17979.00 

576. 1500 

* 1.010000 

23256.00 

626.4400 

1.140000 

28248.40 

648.8400 

■ 1.290000 

33622.93 

657.0667 

■ 1.480000 

39848.54 

648.3400 

1.710000 

46490.44 

625.9360 

m 1.980000 

53950.90 

585.9850 

1 2.260000 

61880.60 

525.7401 

2.820000 

83994.20 

323.7240 

m 3.440000 

111853.9 

131.9580 

■ 4.030000 

138365.6 

2.050000 




Fig. 

A- 7 Sample BASE4 

.TRAJ file 



MACH 

ALTITUDE 

ORBITER 

BOOSTER 

TOTAL 



(FT) 

(LBS) 

(LBS) 

(LBS) 


0.6000000 

9065.000 

1223333. 

2607000. 

6437333. 


0.7400000 

13447.70 

1091033. 

2474000. 

6039033. 


0 . 8700000 

17979.00 

1047400 . 

2369633 . 

5786666 . 


1.010000 

23256.00 

3656867 . 

956333.2 

5569533. 


1 . 140000 

28248.40 

1039534. 

2283134. 

5605801 . 


1 .290000 

33622.93 

1187845. 

2369200. 

5926245. 


1 . 480000 

39848.54 

1296000. 

2461623. 

6219245. 


1 .710000 

46490.44 

1372213. 

2522293. 

6416800. 


1 . 980000 

53950.90 

1385234. 

2554483. 

6494200. 


2.260000 

61880.60 

1395134. 

2539000. 

6473133. 


l . 820000 

83994.20 

1403560. 

2211054. 

5825667. 


- . -v 4 ' 0 0 0 

111853.9 

1408520. 

1810740. 

5030000. 


•. . 0 30000 

1 3836 5.6 

1413240. 

1429796. 

4272833. 



Fig. A 

-8 Sample BASE4. THRUST File 
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LMSC-HEC TR F268592 


MACH 


0.6000000 
0.7400000 
0.8700000 
1.010000 
1. 140000 
1.290000 
1.480000 
1.710000 
1.980000 
2.260000 
2.820000 
3.440000 
4.030000 


ALTITUDE 

(FT) 

9065.000 
13447.70 

17979.00 

23256.00 
28248.40 
33622.93 
39848.54 
46490.44 
53950.90 
61880.60 
83994.20 
111853.9 
1 38365.6 


ORB ANGLE 
(DEG) 

3. 171308 
4.704556 
6.289790 
8.135902 
9.882447 
1 1 .76268 
13.94065 
16.26426 
18.87423 
21.64837 
29.38461 
39. 13108 
48.40594 


Fig. A-9 Sample BASE4. PLUME File 
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RB ANGLE 
(DEG) 

6.072794 

9.008839 

12.04443 

15.57958 

18.92407 
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BASE4 PROGRAM LISTING 


CALULATION OF BASE PRESSURE FOR SSLV w/LRB 


WRITTEN BY: Richard Reheuser 


AT: Lockheed Missiles and Space Co. 

4800 Bradford Drive 
Huntsville, Alabama 38807 

PH: (208)722-4238 


INTEGER INLRB , IN0R8 , IWORB , ITORB , IEL , ILRB 
REAL SREF , MACH # ALT , Q , PI , THTOT , THLRB , THORB , THET , 

PP , PPl , PPI (21) , ALTI (21) , 

DBLRB , DELRB , DC YLRB , XJLRB , DXJLRB , LRB , 

K0LRB , KLRB , EELRB , DJLRB , PCLRB , XSDLRB , 

ANLRB , L 1 LRB , L2LRB , 

DBORB , DEORB , X JORB , DX JORB „ XSORB , TCW , TCT , 

K0ORB , KORB , EEORB , D JORB , PCORB , XSDORB , 

ANORB , L 1 ORB , L20RB , L30RB , 

DBET , DEET , XSET , XSDET , ANET , K0ET , KET , EEET , 

DJET,PCET, 

ABLRB , AELRB , AEFLRB , ACYLRB , AFLRB , PRLRB , PLRBl , 

PLRB2 , CTLRB , TAALRB , CPBLRB , CDBLRB , DFBLRB , 

ABORB , AEORB , AEFORB , PBWORB , PBTORB , CTORB , TA AORB , 

CPBORB , CDBORB , DFBORB , 

ABET , AEFET , CTET , TAAET , CPBET f CDBET , DFBET , 

PBLRB1 , PBLRB2 , PBLRB3 , PBLRB4 , PBLRBA , PBLRB6 , PBLRBT , 

PBORB 1 , PB0RB2 , PB0RB3 , PB0RB4 , PBORBA , PBORBB , PBORBT , 

PBETl , PBET2 , PBET3 , PBET4 , PBETA , PBETB , PBETT , 

DFBLRB1 , DFBORB1 , DFBETl , DFTOTl , PBLRB4A , PB0RB4A , PBET4A , 

DFNLRB , DFNORB , DFNET 
CHARACTER* 1 LINE 

CHARACTER*20 FNAMEl , FNAME2 , FNAME3 , FNAME4 , FNAME6 , FNAME6 
DATA SREF, PI /2690. ,3 . 141692664/ 

DATA ALTI / 0., 10000., 20000., 30000., 40000., 

caa idux assaa oaduxa Qooaa 

DWXrW . , OVvW • , / KVPP. , DvKWu . , WWV. , 

100000 . . 1 10000 . . 1 20000 . . 1 30000 . . 140000 . , 

160000. . 160000. .170000. .180000. .190000. , 

200000./ 

DATA PPI/1 . 10326E+3 , 6 . 96946E+2 , 4 . 86006E+2 , 3 . 01486E+2 , 1 . 88230E+2 , 
1 . 16641E+2 , 7 . 23120E+1 , 4 . 48767E+1 , 2 . 80164E+1 , 1 , 76096E+1 , 
1 . 11428E+1 , 7 . 10413 ,4.89717 ,3.02101 ,2.01393 , 

1 . 36070 , 9 . 29799E-1 , 6 . 36667E-1 , 4 . 39927E-1 , 2 . 94703E-1 , 

1 . 97960E-1/ 

DATA INORB, IWORB, ITORB /3,2,1 / 

DATA PBLRB1 , PBLRB2 , PBLRB3 /0. 0,0. 0,0.0/ 

DATA PBLRB4,PBLRB6, PBLRBT /0. 0,0. 0,0.0/ 

DATA PBORBl , PB0RB2 , PB0RB3 /0. 0,0. 0,0.0/ 

DATA PB0RB4,PB0RB6, PBORBT /0 .0,0 .0,0 .0/ 
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DATA DBORB, DEORB, AFORB /23 . 6 f 7 .5, . 333/ 

DATA K0ORB , EEORB , PCORB /. 0000963 , 7 7 . 5 , 427680 . / 

DATA XJORB, XSORB /9. 2, 10.0/ 

DATA TCW, TCT /.06, .06/ 

C 

DATA PBET1 , PBET2 , PBET3 /0 . 0,0 . 0,0 . 0/ 

DATA PBET 4 , PBET6 , PBETT /0 . 0,0 .0,0 . 0/ 

DATA DBET,XSET, AFET /27 .6,41 . 76 , . 333/ 

C 

DATA FNAME2 / *BASE4 .TRACE*/ 

DATA FNAME3 / » BASE4 . OUT » / 

DATA FNAME4 /> BASE4 . PLUME * / 

DATA F NAMES / , BASE4 . TRAJ •/ 

DATA FNAME6 / * BASE4 . THRUST * / 

C 

WRITE(6, *) 

S WRITE (6, 1000) 

READ (6 , 1001 , ERR=5) FNAME1 
WRITE (6,#) 

C 

OPEN (19 , FILE=FNAME1 , STATUS* * OLD * ) 

OPEN (20 , FILE=FNAME2 , STATUS* ’ NEW > ) 

OPEN (21 , FILE=FNAME3 , STATUS* * NEW* ) 

OPEN (22 , FILE=FNAME4 , STATUS* * NEW * ) 

OPEN (23 , FILE=FNAME6 , STATUS* * NEW * ) 

OPEN (24 , FILE=FNAME6 , STATUS* * NEW * ) 

C 

READ (19, 1001, ERR=998) LINE 
READ (19, 1001 ,ERR*998) LINE 
READ (19, 1001, ERR=998) LINE 
READ (19, 1001 ,ERR=998) LINE 
READ (19, 1001 ,ERR*998) LINE 
READ (19, 1001, ERR=998) LINE 
READ (19, * ,ERR=998) INLRB 
READ (19, *,ERR=998) DBLRB 
READ (19 , * , ERR=998) DELRB 
READ (19, *,ERR=998) DCYLRB 
READ (19,*, ERR=998) XJLR8 
READ (19,*, ERR=998) XSLRB 
READ (19, * ,ERR=998) K0LRB 
READ (19,*, ERR=998) EELRB 
READ (19,*, ERR=998 ) PCLRB 
READ (19,*, ERR=998) K0ET 
READ (19, * , ERR=998) EEET 
READ (19,*, ERR=998) PCET 
READ (19,*, ERR=998) AFLRB 
READ (19, 1001 , ERR*998) LINE 
READ (19, 1001, ERR=998) LINE 
READ (19, 1001, ERR=998) LINE 
READ (19, 1001, ERR=998) LINE 
READ (19, 1001 , ERR=998) LINE 
C 

C LRB CALCULATIONS 

C 

ABLRB = PI * (DBLRB/2 . 0 ) **2 

AELRB = (PI * (DELRB/2 .0) **2) * FLOAT (INLRB) 

AEFLRB = ABLRB - AELRB 
ACYLRB = PI * (DCYLRB/2 .0) **2 
ARLRB = ABLRB/ACYLRB 

L1LRB = XSLRB ♦ DELRB ! NOTE - LI AND L2 ARE SET UP HERE FOR 

L2LRB = XSLRB ♦ DELRB ! A 4 NOZZLE CONFIGURATION, THIS SHOULD 
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ANLRB = L1LRB . 1. 2LRB - AELRB ! BE CHANGED FOR A 2, 3, OR 

XSDLRB = XSLRB/DELRB 


ORB CALCULATIONS 


C 


C 


C 



ABORB = PI • (DBORB/2 . 0) **2 

AEORB = (PI • (DEORB/2.0)* *2) • FLOAT (INORB) 

AEFORB = ABORB - AEORB 

LIORB = XSORB + DEORB 

L20RB = XSORB + DEORB 

L30RB = SQRT (L10RB**2 - (L20RB/2 .0) **2) 

ANORB = 0.6 * L20RB • L30RB - AEORB 
XSDORB = XSORB/DEORB 


ABET = PI * (DBET/2.0)**2 
AEFET = ABET ♦ ABORB ♦ 2* (ABLRB) 
AEET = 3.0 * AELRB 

DEET = 2.0 * SQRT (AEET/PI) 

ANET = 0.6 * ABET 

XSDET = XSET/DEET 


WRITE(20,*) 
WRITE (20, 2006) 
WRITE (20, 2007) 
WRITE (20, 2008) 
WRITE (20, 2009) 
WRITE (20, 20X0) 
WRITE(20,2011) 
WRITE (20, 2012) 
WRITE (20, 2013) 
WRITE (20, 2014) 


ABLRB, ABORB, ABET 
AELRB, AEORB, AEET 
AEFLRB, AEFORB, AEFET 
ACYLRB, LIORB, DEET 
ARLRB , L20RB , ANET 
L1LRB.L30RB, XSDET 
L2LRB, ANORB 
ANLRB, XSDORB 
XSDLRB 


WRITE (20,*) 
WRITE(20,*) ’**■ 




WRITE(20,*) **' 
WRITE (20,*) 
WRITE (21, 2000) 
WRITE (21, 2001) 
WRITE(21 ,*) 
WRITE (21, 2002) 
WRITE (21,*) 
WRITE (21, 2003) 
WRITE (21, 2004) 



ORB ANGLE 

WRITE (22,*) ’ 

» 

MACH 

RB ANGLE * 

ALTITUDE 

WRITE (22,*) ’ 

» 

WRITE (22,.) 

(DEG) » 

(FT) 

(DEG) 

WRITE (23,*) ’ 

MACH 

ALTITUDE 

Q* 

WRITE (23,*) ’ 
WRITE (23,*) 


(FT) 

(PSF) * 

WRITE (24, ») * 

MACH 

ALTITUDE 

ORBITER 

* BOOSTER 

TOTAL’ 


WRITE (24,.) • 


(FT) 

(LBS) 

> 

WRITE (24 , *) 

(LBS) 

(LBS) ’ 



BEGIN TRAJECTORY LOOP 


READ (1 9 , * , ERR=998 , END=999 , IOSTAT=IER) MACH , ALT , Q , THTOT , THLRB 
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THORB = THTOT - 2 0*THI RB 

THET = THTOT 

CTLRB = THLRB/ (QaABLRB) 

CTORB = THORB/ (Q*ABORB) 

CTET = THET/ (Q.AEFET) 

TAALRB = THLRB/ABLRB 
TAAORB = THORB/ABORB 
TAAET = THET/AEFET 

UNCORRECTED BASE PRESSURE RATIO - (ALL M; ALL ELEMENTS) 

IF (MACH. LT. 1.8) THEN 
CALL EXTCT (MACH, CTLRB, PBLRB1A) 

CALL EXTCT(MACH, CTORB, PBORB1A) 

CALL EXTCT (MACH, CTET, PBET1A) 

PBLRB1 = PBLRB1A 
PBORB1 = PBORB1A 
PBET1 = PBET1A 
END IF 

IF (MACH. GE. 1.8. AND. MACH. LE. 2. 2) THEN 
CALL EXTCT (MACH, CTLRB, PBLRB1 A) 

CALL EXTCT (MACH, CTORB, PBORB1A) 

CALL EXTCT (MACH, CTET, PBET1 A) 

CALL EXTCA(ALT, TAALRB, PBLRB1B) 

CALL EXTCA(ALT, TAAORB, PBORB1B) 

CALL EXTCA(ALT, TAAET, PBET1B) 

PBLRB1 = (PBLRB1A ♦ PBLRBlB)/2.0 
PBORB1 = (PBORB1A ♦ PBORBlB)/2.0 
PBET1 = (PBET1A * PBET1EI )/2.0 
END IF 

IF (MACH. GT. 2. 2) THEN 
CALL EXTCA(ALT, TAALRB, PBLRB1B) 

CALL EXTCA(ALT, TAAORB, PBOFIB1B) 

CALL EXTCA (ALT, TAAET, PBET1B) 

PBLRB1 = PBLRB1B 
PBORB1 = PBORB1B 
PBET1 = PBET1B 
END IF 

CORRECTION FOR NOZZLE EXTENSION - (ALL M; LRB.ORB) 

DXJLRB = ABS (X JLRB/DBLRB - 0 34) 

CALL EXTDXJ (MACH, DXJLRB, PBLRB2) 

IF (X JLRB/DBLRB. LT. 0.34) PBLRB2 = -PBLRB2 
DXJORB = ABS (X JORB/DBORB - 0.34) 

CALL EXTDXJ(MACH, DXJORB, PB0RB2) 

IF (X JORB/DBORB. LT. 0.34) PB0RB2 = -PBORB2 

CORRECTION FOR NON-CYLINDRICAL - (1.60 .GE. M .LE. 3.6; LRB) 

IF(MACH.LT. 1 .0) THEN 
PLRB1 =0.0 
PLRB2 =0.0 
PBLRB3 =0.0 
ENDIF 

IF (MACH . GE . 1 . 0 . AND . MACH . LT . 1 . 6) THEN 
CALL EXTCYL (ARLRB.MACH.PRLRB) 

PLRB1 = PBLRB1 ♦ PBLRB2 ♦ PBLRB3 
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PLRB2 = PLRB1/PRLRB 

PBLRB3 - PLRB2 - PLRB1 

PBLRB3 = PBLRB3 • ((2.0 * MACH) - 2.0) 

END IF 

IF (MACH . GE . 1 . 6 . AND .MACH . LT . 3 . 6) THEN 
CALL EXTC YL ( ARLRB , MACH , PRLRB) 

PLRB1 = PBLR81 ♦ PBLRB2 ♦ PBLRB3 
PLRB2 = PLRB1 /PRLRB 
PBLRB3 = PLRB2 - PLRB1 
END IF 

IF (MACH. GE. 3. B. AND.MACH.LT. 4.0) THEN 
CALL EXTC YL( ARLRB, MACH, PRLRB) 

PLRBl = PBLRB1 ♦ PBLRB2 ♦ PBLRB3 
PLRB2 = PLRBl /PRLRB 
PBLRB3 = PLRB2 - PLRBl 
PBLRB3 = PBLRB3 • ((-2.0 • MACH) ♦ 8.0) 
END IF 

IF (MACH. GE. 4.0) THEN 
PLRBl = 0.0 
PLRB2 =0.0 
PBLRB3 =0.0 
ENDIF 

CORRECTIONS FOR FINS - (ALL M; ORB) 

CALL EXTF IN (MACH , PBWORB) 

CALL EXTF IN (MACH, PBTORB) 

PBWORB = PBWORB • (TCW/0.1) • (IWORB/4.0) 
PBTORB = PBTORB * (TCT/0.1) * (ITORB/4.0) 
PBORBB = PBWORB ♦ PBTORB 

TOTAL ALL Pb/Pi 

PBLRBT = PBLRB1 ♦ PBLRB2 ♦ PBLRB3 ♦ PBLRB6 
PBORBT = PBORB1 ♦ PB0RB2 ♦ PB0RB3 ♦ PBORBB 
PBETT = PBET1 + PBET2 ♦ PBET3 ♦ PBET6 

BASE PRESSURE COEFF 

CPBLRB = (PBLRBT - 1.0)/(0.7 • MACH*«2) 

CPBORB = (PBORBT - 1.0)/(0.7 * MACH*»2) 

CPBET = (PBETT - 1.0)/(0.7 * MACH**2) 

BASE DRAG COEFF 


C 


C 


CDBLRB = -(CPBLRB • AEFLRB/SREF) 
CDBORB = -(CPBORB • AEFORB/SREF) 
CDBET = -(CPBET • ABET/SREF) 


CORRECTION FOR NOZZLE SPACING - (ALL M; LRB,ORB,ET) 


DJLRB = K0LRB • 
DJORB = K0ORB • 
DJET = K0ET * 


(ALT/1000.) • 
(ALT/1000.) • 
(ALT/1000.) * 


(PCLRB) *«0 . 8 • 
(PCORB) **0 . 8 * 
(PCET) **0 . 8 * 


(EELRB) •* (-0. B) 
(EEORB) ** (-0 . 6) 
(EEET) ** (-0 . B) 


WRITE (22 , *) MACH, ALT, DJORB, DJLRB 

WRITE(23, *) MACH, ALT, Q 

WRITE (24 , *) MACH, ALT, THORB , THLRB , THTOT 


IF (XSDLRB . GT .0 . 0) THEN 
CALL EXTXS1 (DJLRB, XSDLRB, PBLRB4) 
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ENOIF 

C 

IF (XSDORB . GT . 0 . 0) THEN 
CALL EXTXS2 (DJORB, XSDORB ,PB0RB4) 
END IF 

IF (XSDET . GT .0 .0) THEN 
CALL EXTXS2(D JET, XSDET, P0ET4) 
ENDIF 

NOZZLE SPACING EFFECT DRAG 

CALL ATP (ALT , ALTI ,PPI , 21 , PP) 

PP1 = PP * 2.089 

PBLRB4A = PBLRB4 * PPl 
PB0RB4A = PB0RB4 * PPl 
PBET4A = PBET4 * PPl 

DFNLRB = -PBLRB4A * AFLRB * ANLRB 

DFNORB = -PB0RB4A * AFORB * ANORB 

DFNET = -PBET4A * AFET * ANET 


BASE DRAG 


C 


C 


DFBLRB1 = 2.0 * CDBLRB*Q*SREF 

DFBORB1 = CDBORB*Q*SREF 

DFBET1 = CDBET*Q*SREF 

DFBT0T1 = DFBLRB1 «• DFBORB1 ♦ DFBET1 

DFBLRB = DFBLRBl ♦ 2.0 *DFNLRB 

DFBORB = DFBORB1 ♦ DFNORB 

DFBET = DFBETl ♦ DFNET 

DFBTOT = DFBLRB ♦ DFBORB ♦ DFBET 


WRITE (20,*) 
WRITE(20, *) » 
WRITE (20,*) ’ 
WRITE (20,*) * 
WRITE(20 , •) 1 
WRITE (20 , *) * 
WRITE(20 , *) * 
WRITE(20 , *) * 
IF (MACH. LT. 2 
WRITE (20, *) 
WRITE (20, *) 
WRITE (20, *) 
ENDIF 

IF (MACH. GE. 2 
WRITE (20, *) 
WRITE (20, *) 
WRITE (20 , ♦) 
ENDIF 

WRITE (20 , *) ' 
WRITE(20 , *) 
WRITE (20 , *) * 
WRITE (20 , *) * 
WRITE (20 , *) * 
WRITE (20 , *) * 
WRITE (20 , *) ’ 
WRITE (20 , *) 


MACH 

ALT 

Q 

PRESS 

THRUST TOT 
BOOSTER THRUST 
ORB THRUST 


(FT) 

(PSF) 

(PSF) 

(LB) 

(LB) 

(LB) 


* , MACH 

* , ALT 

\Q 

* ,PPi 

► ,THTOT 
* ,THLRB 
* , THORB 


0) THEN 
CTBOOSTER 
CTORB 
CTET 

* , CTLRB 
* , CTORB 
1 , CTET 


0) THEN 
T/A BOOSTER 
’ T/A ORB 
’ T/A ET 

» , TAALRB 
* , TAAORB 
* , TAAET 


PLUME ANGLE 

(RB ORB) 

’ ,DJLRB, DJORB 

UNCORRECTED 
NOZZLE EXT 
NON CYLINDER 
FINS 

TOTAL Pb/Pi 

(RB ORB ET) 
(RB ORB ET) 
(RB ORB ET) 
(RB ORB ET) 
(RB ORB ET) 

» ,PBLRB1,PB0RB1,PBET1 
* ,PBLRB2,PB0RB2,PBET2 
’ , PBLRB3 , PB0RB3 , PBET3 
• , PBLRB5 , PBORBS , PBETS 
* , PBLRBT , PBORBT, PBETT 
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c 

c 

1000 

1001 

1002 

2000 

2001 

2002 

2003 

2004 
2006 

2006 

2007 * 

2008 

2009 

2010 
2011 

2012 

2013 

2014 

C 

998 

999 

C 

C****» 

c 

c 

c 

c 


c 


WRITE (20, *) 
WRITE (20, *) 
WRITE (20,*) 
WRITE (20, *) 
WRITE (20, *) 
WRITE(20, *) 
WRITE(20, *) 
WRITE (20, •) 
WRITE(20, *) 
WRITE(20, *) 
WRITE(20, *) 
WRITE(20, *) 
WRITE(20, *) 
WRITE(20, *) 


' CPB 


(RB ORB 

* CDB 


(RB ORB 

’ DFB 

(wo/NS) 

(RB ORB 

* DFB 

(wo/NS) 

(TOTAL) 

* NOZ 

SP Pb/Pi 

(RB ORB 

’ NOZ 

SP Pb 

(RB ORB 

' NOZ 

SP DRAG 

(RB ORB 

* DFB 

(w/NS) 

(RB ORB 

> DFB 

(w/NS) 

(TOTAL) 


ET) » , CPBLRB , CPBORB , CPBET 
ET) * , CDBLRB , CDBORB , CDBET 
ET) > , DFBLRB1 , DFBORBl ,DFBET1 
, ,DFBTOTl 

ET) * , PBLRB4 , PB0RB4 , PBET4 
ET) * , PBLRB4A, PBORB4A , PBET4A 
ET) > , DFNLRB , DFNORB , DFNET 

ET) * , DFBLRB , DFBORB , DFBET 
* ,DFBT0T 


WRITE (2 1 , 2006 ) MACH , ALT , Q , DFBTOT , DFBORB , DFBLRB , DFBET 
GO TO 100 


FORMAT (’ ENTER INPUT FILE NAME — > \i) 

FORMAT (A) 

FORMAT ( * ENTER BOOSTER TYPE »,/, 

• 1-SRB 2-MLRB1 3-MLRB2 4-GDLRB3 6-GDLRB4 — > »,$) 
FORMAT (2 (/)) 

FORMAT (26X , 8 ( * X * ) , IX , 8 ( * X * ) ,2X, * (XXX) ») 

FORMAT ( 30X , ’ TOTAL BASE DRAG » ) 

FORMAT (6X, *MACH ALT (FT) Q(psf) LV(lb*) \ 

* ORB(lbs) 2RB ( I bs) ET(lbs)*) 

FORMAT (3X , 74 ( * - * ) ) 

FORMAT (6X , F4 . 2 , 2X , F10 . 2 , 2X ,F6 . 1 , 2X, F10 . 2 , 2X ,F10 . 2 , 2X , 

F10 . 2 # 2X, F10 . 2) 

FORMAT (T2 , 1 ABLRB * , T13,F11 . 5 ,T28, * ABORB \T38,FU.S, 

T63 , ’ABET \T63,F11.5) 

FORMAT (T2 , *AELRB * ,T13,F11 . 5,T28, ’AEORB \T38,F11.5, 

T53, * AEET \T63,F11.6) 

FORMAT (T2 , »AEFLRB\T13 ,F11 . 5,T28, *AEFORB* ,T38,F11 .6, 

T63 , *AEFET \T63,F11.6) 

FORMAT (T2 , ’ ACYLRB* ,T13,Fll .5,T28, ’LIORB \T38,F11.S, 

T53, ’DEET \T63,F11.6) 

FORMAT (T2 , * ARLRB * ,T13,F11 .5,T28, »L20RB \T38,F11.5, 

T53, *ANET \T63,F11.6) 

FORMAT (T2 , *L1LRB * ,T13 ( F11 .6,T28, ^SORB \T38,F11.6, 
T63/XSDET \T63,F11.6) 

FORMAT (T2 , 'L2LRB > , T13 , Fll . 6, T28 , » ANORB >,T38,F11.6) 
FORMAT (T2 , 'ANLRB * ,T13,FU .5,T28, ’XSDORB* ,T38,Fll .5) 
FORMAT (T2 , ’XSDLRB* ,T13,F11 .6) 

WRITE (6 , *) ALT, MACH, Q, THTOT, THLRB 
WRITE(6, *) * READ ERROR * , IER 
END 


SUBROUTINE EXTCT (MACH , CTHRUST , PBB7) 

EXTRAPOLATE Pb/Pi vs M and CT PB(M,CT) 

REAL M(9) , CT (7) ,PB(9,7) , MACH , CTHRUST, PERM, PERCT , 

PBB1 , PBB2 , PBB3 , PBB4 , PBB6 , PBB6 , PBB7 

DATA M /0. 00 ,0.60, 0.90,1 .00,1.10,1.25,1 .60,1 .75,2.00/ 
DATA CT / 0.0, 5.0,10.0,16.0,20.0,25.0,30.0/ 
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c 


c 

20 


C 

30 


C 


c 

c 

c 

c 

c 

c 


c 


c 


DATA PB /0. 9600, 0.9500, 0.9250, 0.8800, 0.8600, 0.8000, 0.7400, 
0.6000,0.4600, 

0.9600,0.9600,0.9260,0.8900,0.8800,0.8750,0.9100, 

0.8900,0.8700, 

0 . 9500 , 0 . 9600 ,0.9250,0. 9000 , 0 . 8950 ,0.9250,1.0400, 
1.1000,1.1600, 

0.9500,0.9500,0.9250,0.9150,0.9100,0.9500,1 .0700, 
1.1900,1.3100, 

0.9500,0.9600,0.9250,0.9260,0.9200,0.9700,1.1000, 

1.2400,1.3800, 

0.9600,0.9600,0.9260,0.9400,0.9400,0.9850,1.1300, 

1.2900,1.4600, 

0 . 9500 , 0 . 9600 ,0.9260,0. 9600 , 0 . 9500 ,0.9860,1.1 500 , 
1.3200,1.4900/ 

DO 1=2,9 

IF (M (I) . GT . MACH) THEN 
IM = I 
GO TO 20 
ENDIF 

IF (I . EQ . 9) IM = 9 
ENDDO 

DO J=2 , 7 

IF(CT(J) .GT.CTHRUS T) THEN 
ICT = J 
GO TO 30 
ENDIF 

IF ( J . EQ . 7) ICT = J 
ENDDO 

PERM = (M (IM) -MACH) / (M(IM) -M(IM-l) ) 

PERCT = (CT (ICT) -CTHRUST) / (CT (ICT) -CT (ICT-1) ) 

PBBl = PB(IM-1, ICT-1) 

PBB2 = PB(IM, ICT-1) 

PBB3 = PB (IM-1 , ICT) 

PBB4 = PB(IM, ICT) 

PBB6 = PBB2 - (PBB2 - PBB1)*PERM 

PBB6 = PBB4 - (PBB4 - PBB3) *PERM 

PBB7 = PBB6 - (PBB6 - PBB5) *PERCT 

RETURN 

END 


SUBROUTINE EXTCA (ALT , TAAA , PB07) 

EXTRAPOLATE Pb/Pi vs C/A and ALT PB (TAA , H) 

REAL TAA (6) ,H(14) ,PB(6,14) , ALT, TAAA, PERTAA,PERH, 

PBBl , PBB2 , PBB3 , PBB4 , PBB5 , PBB6 , PBB7 

. DATA TAA / 3000 ., 3600, 4000 5000 . ,6000 ., 10000 . / 

DATA H / 0.0, 10000.0, 20000.0, 30000.0, 40000.0, 50000.0 

60000.0, 70000.0, 80000.0, 90000 .0,1 00000 .0,11 0000 0 * 
120000.0,13 0000.0/ ' ' 

DATA PB /0 50,0.65,0.60,0.73,0.78,0.88, 
0.58,0.65,0.70,0.80,0.90,0.98, 
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c 


c 

20 


C 

30 


C 


C 

C. 

C 

C 

c 

c 


c 


c 


0.63, 0.70, 0.78, 0.90,1 .00,1 .08, 
0.70,0.76,0.86,1.00,1.10,1.20, 
0.78,0.86,0.93,1.11,1.23,1.30, 
0.85,0.93,1.00,1.20,1.33,1.40, 
0.90,0.98,1.13,1.33,1.45,1.50, 
0.98,1.13,1.30,1.55,1.66,1.73, 
1.08,1.28,1.50,1.80,1.90,1.95, 
1.18,1.43,1.70,2.03,2.25,2.20, 
1.30,1.60,1.93,2.30,2.46,2.48, 
1.43,1.80,2.18,2.60,2.75,2.85, 
1.65,2.00,2.45,2.95,3.20,3.36, 
1.68,2.30,2.85,3.45,3.85,4.10/ 

DO 1=2,6 

IF(TAA(I) .GT.TAAA) THEN 
ICA = I 
GO TO 20 
END IF 

IF(I.EQ.6) ICA = I 
ENDDO 

DO J=2 , 14 

IF (H(J) .GT.ALT) THEN 
IH = J 
GO TO 30 
ENDIF 

IF ( J . EQ. 14) IH = J 
ENDDO 

PERTAA = (TAA (ICA) -TAAA) / (TAA (ICA) -TAA (ICA-1) ) 
PERM = (H (IH) -ALT) / (H (IH) -H (IH-1) ) 

PBB1 = PB (ICA-1 , IH-1) 

PBB2 = PB (ICA, IH-1) 

PBB3 = PB (ICA-1 , IH) 

PBB4 = PB (ICA, IH) 

PBB6 = PBB2 - (PBB2 - PBB1) *PERTAA 

PBB6 = PBB4 - (PBB4 - PBB3) sPERTAA 

PBB7 = PBB6 - (PBB6 - PBB5)*PERH 

RETURN 

END 


SUBROUTINE EXTDXJ (MACH , DELXJ , DPBB7) 

EXTRAPOLATE DPb/Pi vs M and DXJ DPB(M,DXJ) 

REAL M(18) ,DXJ(11) ,DPB(18,11) , MACH , DELXJ , PERM, PERXJ , 

DPBBl , DPBB2 , DPBB3 , DPBB4 , DPBB5 , DPBB6 , DPBB7 

DATA M /0. 0,0. 4, 0.8, 1.2, 1.6, 2. 0,2. 4, 2. 8, 3. 2, 3. 6, 4. 0,4. 4, 
4. 8, 5. 2, 5. 6, 6. 0,6. 4, 6. 8/ 

DATA DXJ / 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0/ 

DATA DPB /0. 000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 

0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 0 • 000 . 0 • 000 , 
0.000,0.000,0.000,0.000, 

0.000,0.000,0.000,0.007,0.024,0.032,0.028, 

0 . 020 ,0.012,0. 007 ,0.003,0. 000 , 0 . 000 , 0 . 000 , 

0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 
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C 

20 


C 

30 


C 


c 

c 

c 

c 


0.000,0.005,0.006,0.016,0.048,0.062,0.066, 
0.046,0.032,0.024,0.016,0.008,0.002,0.000, 
0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 
0.010,0.012,0.014,0.032,0.080,0.094,0.086, 
0.072,0.066,0.040,0.030,0.018,0.006,0.000, 
0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 
0.020,0.022,0.024,0.063,0.116,0.132,0.120, 
0.100, 0.080, 0.063, 0.046, 0.028, 0.012, 0.003, 
0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 
0.029,0.031,0.038,0.077,0.166,0.168,0.160, 
0.128,0.104,0.083,0.062,0.040,0.022,0.008, 
0 . 000 , 0 . 000 , 0 . 000 , 0 . 000 , 
0.039,0.040,0.062,0.103,0.188,0.204,0.184, 
0.166,0.129,0.104,0.078,0.064,0.034,0.015, 
0 . 002 , 0 . 000 , 0 . 000 , 0 . 000 , 
0.047,0.060,0.066,0.126,0.224,0.240,0.216, 
0.184,0.163,0.126,0.096,0.068,0.044,0.024, 
0.009,0.000,0.000,0.000, 
0.066,0.060,0.080,0.163,0.264,0.276,0.244, 
0.212,0.179,0.147,0.114,0.084,0.066,0.032, 
0.016,0.003,0.000,0.000, 
0.066,0.070,0.102,0.184,0.300,0.316,0.276, 
0.240,0.203,0.161,0.130,0.098,0.066,0.041, 
0.023,0.007,0.000,0.000, 
0.077,0.081,0.130,0.213,0.336,0.348,0.308, 
0.268,0.229,0.187,0.146,0.110,0.076,0.049, 
0.030,0.013,0.000,0.000/ 


DO 1*2 , 18 

IF(M(I) . GT.MACH) THEN 
IM : I 
GO TO 20 
ENDIF 

IF (I . EQ. 18) IM = 18 
ENDDO 

DO J*2,ll 

IF(DXJ(J) .GT.DELXJ) THEN 
IXJ = J 
GO TO 30 
ENDIF 

IF(J.EQ.ll) IXJ s J 
ENDDO 

PERM = (M (IM) -MACH) /(M(IM)-M(IM-l)) 

PERXJ = (DX J(IXJ)-DELXJ)/(DXJ(IXJ) -DXJ (IXJ-1) ) 
DPBB1 = DPB(IM-1, IXJ-1) 

DPBB2 = DPB(IM, IXJ-1) 

DPBB3 = DPB(IM-1,IXJ) 

DPBB4 = DPB(IM.IXJ) 

DPBB6 = DPBB2 - (DPBB2 - DPBB1) (PERM 

DPBB6 = DPBB4 - (DPBB4 - DPBB3) *PERM 

DPBB7 = DPBB6 - (DPBB6 - DPBB6)*PERXJ 

RETURN 

END 


SUBROUTINE EXTCYL (AR, MACH ,PR) 
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c 

C 

c 

c 


c 


c 

20 


C 


C 

C 

C 

c 

c 

c 

c 

c 


c 


c 


c 


c 

20 


c 

30 


EXTRAPOLATE Pressure Ratio of NON-CYL vs 
Mach and Base Area Ratio PR = f(M,AR) 

REAL M(7) ,S(7) ,AR,MACH,PR 

DATA M /0. 00, 1.60, 1.76, 2. 00, 2. 50, 3. 60, 30. 00/ 
DATA S /. 000, .000, .060, .140, .105, .080, .080/ 

DO 1=2,7 

IF (M(I) . GT . MACH) THEN 
IM = I 
GO TO 20 
ENDIF 

IF (I .EQ . 7) IM = 7 
ENDDO 

PERM = (M (IM) -MACH) / (M (IM) - M(IM-l) ) 

SLOPE = S (IM) - (S (IM) - S (IM-1) ) *PERM 
PR = SLOPE • (AR - 1.0) ♦ 1.0 

RETURN 

ENO 


SUBROUTINE EXTXS1 (D J, XSDD , PBB7) 

EXTRAPOLATE Pb/Pi vs DJ and XSD PB(PL,XSD) 
********* FOUR NOZZLE CONFIGURATION ****** 


REAL PL (10) ,X$D (4) ,PB(10,4) ,DJ,XSDD,PERPL,PERXSD, 
PBBl , PBB2 , PBB3 , PBB4 , PBB5 , PBB6 , PBB7 


DATA PL /0. 0,10. 0,20. 0,30. 0,40. 0,50. 0,60. 0,70. 0,80. 0,90.0/ 
DATA XSD/1. 00,1. 22, 1.36, 1.48/ 


DATA PB / 0.0, 


0. 

.3, 

2.0, 

7, 

•0, 

12.0 

0. 

0, 

0.0, 

1 

■0, 

4.8 

0, 

■0, 

0.0, 

0 

•0, 

1.0 

0 

>0, 

0.0, 

0 

0, 

0.2 


9.6.14.6.19.5.24.6.29.6, 
5.0,10.3,15.3,20.3,25.3, 

1.6, 6.5,10.5,15.5,20.6/ 


DO 1=2,10 

IF (PL (I) . GT.D J) THEN 
IPL = I 
GO TO 20 
ENDIF 

IF (I . EQ . 10) IPL = 10 
ENDDO 


DO J=2 , 4 

IF (XSD (J) .GT.XSDD) THEN 
IXSD = J 
GO TO 30 
ENDIF 

IF ( J . EQ . 4) IXSD = J 
ENDDO 


PERPL = (PL (IPL) “D J) / (PL (IPL) -PL (IPL-1)) 

PERXSD = (XSD (IXSD) -XSDD) / (XSD (IXSD) -XSD ( IXSD- 1) ) 
PBBl = PB (IPL-1 , IXSD-1) 
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PBB2 = PB (IPL , IXSD-1) 

PBB3 = PB(IPL-1 , IXSD) 

PBB4 = PB (IPL , IXSD) 

PBB6 = PBB2 - (PBB2 - PBB1)*PERPL 

PBB6 = PBB4 - (PBB4 - PBB3)*PERPL 

PBB7 = PBB6 - (PBB6 - PBB6)*PERXSD 

IF (PBB7 . LT . 0 . 0) PBB7 =0.0 
IF (PBB7 . GT . 8 . 0) PBB7 = 8.0 

RETURN 

END 


SUBROUTINE EXTXS2 (D J , XSDD , PBB7) 

EXTRAPOLATE Pb/Pi vs DJ and XSD PB(PL,XSD) 

********* THREE NOZZLE CONFIGURATION ********* 

REAL PL (10) , XSD (3) ,PB(10,3) ,DJ,XSDD, PERPL , PERXSD , 

PBBl , PBB2 , PBB3 , PBB4 , PBB5 , PBB6 , PBB7 

DATA PL /0. 0,10. 0,20. 0,30. 0,40. 0,60. 0,60. 0,70. 0,80. 0,90.0/ 
DATA XSD/1 .00,1 .22,3.00/ 


DATA PB / 0.0, 

0.0, 

0.0, 

0.0, 0.3, 1.1, 

3.6, 

8.6,13.6,18.6, 

0.0, 

0.0, 

0.0, 

0.0, 0.1, 0.6, 

1.8, 

6.0,11.0,16.0, 

0.0, 

0.0, 

0.0, 

0.0, 0.0, 0.1, 

0.4, 

3.0, 7.8,12.8/ 


DO 1=2,10 

IF (PL (I) .GT.DJ) THEN 
IPL = I 
GO TO 20 
END IF 

IF (I . EQ . 10) IPL = 10 
ENDDO 
C 

20 DO J=2 , 3 

IF (XSD (J) .GT.XSDD) THEN 
IXSD = J 
GO TO 30 
ENDIF 

IF ( J . EQ .3) IXSD = J 
ENDDO 
C 

30 PERPL = (PL (IPL) -D J) / (PL (IPL) -PL (IPL-1) ) 

PERXSD = (XSD (IXSD) -XSDD) / (XSD (IXSD) -XSD (IXSD -1)) 
PBBl = PB (IPL-1 , IXSD-1) 

PBB2 = PB(IPL, IXSD-1) 

PBB3 = PB (IPL-1, IXSD) 

PBB4 = PB (IPL, IXSD) 

PBB6 = PBB2 - (PBB2 - PBB1)*PERPL 

PBB6 = PBB4 - (PBB4 - PBB3)*PERPL 

PBB7 = PBB6 - (PBB6 - PBB5) *PERXSD 

IF(PBB7 .LT.0.0) PBB7 = 0.0 
IF (PBB7 . GT .8.0) PBB7 = 8.0 
C 

RETURN 

END 
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c 

SUBROUTINE EXTFIN (MACH , PBB3) 

EXTRAPOLATE Pb/Pi for Fins vs M PB(M) 

REAL M (9) , PB (9) , MACH , PBBl ,PBB2 , PBB3 

DATA M / 0.00,0.60,0.76,1.00,1.26,1.60,2.00,3.00,30,00/ 
DATA PB / 0 . 000 , -0 . 006 ,-0.012,-0. 030 ,-0.044,-0.049, 
-0.050, >0.048, -0.048/ 


DO 1=2,9 

IF (M (I) . GT . MACH) THEN 
IM = I 
GO TO 20 
ENDIF 

IF (I . EQ . 9) IM = 9 
ENDDO 

PERM = (M (IM) -MACH) / (M (IM) -M(IM-l) ) 
PBBl = PB(IM-l) 

PBB2 = PB (IM) 

PBB3 = PBB2 - (PBB2 - PBB1)*PERM 

RETURN 

END 


SUBROUTINE ATP (ALT , ALTI , PPI ,N, PP) 
DIMENSION ALTI (21) , PPI (21) 

DO 10 1=2, N 

IF (ALT. LE. ALTI (I)) GO TO 20 
10 CONTINUE 
I=N 

20 PCT= (ALT-ALTI (1-1) ) / (ALTI (I) -ALTI (1-1) ) 
PP=PPI (1-1) +PCT* (PPI (I) -PPI (1-1 ) ) 

RETURN 

END 
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